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1. Introduction 

A glass may simply be defined as a non-crystalline solid and is described by the 
predominance of local order. Translational symmetry, i. e. long range order, is lacking. 
The term frozen liquid stresses the structural similarity to a liquid. 

By glass transition we define the process of solidification of a supercooled liquid to 
a non-crystalline solid, more exactly we may call it a liquid-to-glass transition. When 
crystallization during cooling can be avoided the structural a-relaxation time r increases 
dramatically. In a normal liquid the relaxation time is of the order of 1 psec, it increases 
continuously within a narrow temperature range over many orders of magnitude to 
attain values in the range of seconds, hours or even days. Then the supercooled liquid 
is conventionally called a glass. 

Many theoretical concepts have been proposed to explain the glass transition among 
these free-volume and thermodynamic theories. In the free-volume theories, introduced 
by Cohen and Turnbull ^ and substantiated by Cohen and Crest [0 the slowing 
down of the dynamics is related to the decrease of free volume and excess entropy. 
Thermodynamic theories first formulated by Cibbs and Di-Marzio and Adams and 
Gibbs m assume an underlying second-order phase transition in the limit r ^ oo to 
account for the Kauzmann paradox. 

A new approach to the problem was opened up after the mode-coupling 
approximation has been introduced first by Leutheusser [^] and independently by 
Bengtzelius, Cotze and Sjogren [0. It provides an appproximation scheme which allows 
for the calculation of the evolution of time correlations from a microscipic equation of 
motion and basically deals with the two well known phenomena: cage and backflow 
effect. In order to understand the structural relaxation it is proposed to investigate the 
short-time limit r ~ 1 psec where relaxational motion begins to evolve from elementary 
vibrations of atoms or molecules. 

Even in its most reduced schematic form mode-coupling theory (MCT) is able to 
describe essential features of structural relaxation and provides new predictions which 
can be tested experimentally. Numeric solutions of mode-coupling equations describe 
the dynamics of computer liquids or colloids with an accuracy of better than 10% without 
any adjustable parameters [0, ||, ^. The structure of the theory suggests that the 
asymptotic results hold also for more complicated "real" glass fromers. However, only 
experiments can decide to which extent "real" systems reach the asymptotic behaviour 
predicted by theory. An increasing number of measurements on many different materials 
has been devoted to this question and give evidence for the relevance of the mode- 
coupling theory |^ . 



In the present review, we survey results for one often studied model system, the 
molecular liquid orthoterphenyl (OTP) (section]^). The glass former orthoterphenyl can 
be considered to be a very good candidate to test the validity of such a theory as the 
molecules are relatively symmetric and rigid. The model character is underlined by the 



recurrent use for experimental studies and tests of various theories 
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Neutron scattering is particularly well suited to test microscopic theories as it 
measures directly the density and tagged-particle correlations on time scales of psec 
. . . nsec and on length scales of intermolecular distances and therefore provides an 
ideal link between theory and experiment. Thus, neutron scattering experiments have 
been performed in a large variety of different glass forming systems like ionic glasses 
|T^, |T^, |T^, |T^, polymers ||TB|, |rH| and even metallic glasses [|TU], ^ A comparison 
can be found in |^ . 

To establish the context, we first describe some glass transition signatures 
(section ^ and some examples of the evolution of structural relaxation on time 
scales T ^ Insec accessible by conventional techniques. After some remarks on the 
experimental method neutron scattering (section ^) a short introduction into the mode- 
coupling theory of the glass transtion (section ^ is given. We then concentrate on the 
fast dynamics on the psec. . .nsec time scale. Tagged-particle and collective motion in 
liquid, glassy and crystalline (section |) orthoterphenyl have been investigated as a 
function of temperature (section H), pressure (section |^) and scattering vector. In this 
range time correlation functions evolve from vibrations (section ^ at very short times 
to structural relaxation through an intermediate regime for which the mode-coupling 
theory makes the most prominent predictions. 

The neutron scattering experiments on OTP are probably the most intensive study 
on a molecular glass former up to now. Here, we not only summarize previously 
published results, but we also present many new yet unpublished data. Some 
experimental limitations are discussed in detail. Those figures, which are here published 
for the first time, or where we include new unpublished data are denoted by "*" in the 
figure caption. All other figures (or data) are taken from the literature as stated in the 
caption. This review marks the end of a long term project that was supported over more 
than 10 years by BMBF under project numbers 03si2mai, 03fu4dor4 and 03fu5dor1. 
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2. Orthoterphenyl 

Every liquid may be solidified as a glass provided crystallization can be bypassed. 
Practically, atomic liquids cannot be supercooled into the region of interest without 
immediate crystallization. The discovery of a concentration driven glass transition in 
colloidal suspensions provided new model systems in which the dynamics on the length 
scale of the nearest neighbours is probed by dynamic light scattering. Indeed, very 
impressive confirmations of MCT predictions have been provided by light scattering 
studies in such systems 0, ||, H. 

According to mode-coupling theory the dramatic increase of the relaxation time 
with decreasing temperature or increasing pressure is due to relatively subtle but smooth 
changes in the static structure factor. Since the theory was originally derived for "simple 
liquids", it can be expected to be more reliable for simple systems without internal 
degrees of freedom than for more complicated molecules. 

Molecular liquids form stable supercooled liquids in favorable cases. However, 
processes related to the molecular complexity which allows for supercooling should not 
screen phenomena that are predicted by MCT, which, originally, has been developed for 
simple atomic liquids with hard-core or Lennard- Jones interaction potentials without 
internal degrees of freedom. 

The van der Waals liquid Orthoterphenyl (OTP: 1,2-diphenylbenzene, C18H14, 
boiling point Tb = 605 K, melting point Tin=329 K) is extensively studied for more 
than 40 years as an archetypical glass former by many methods and is representative 
for a whole class of glasses named fragile EH], whose characteristic is an extremely 



non-Arrhenius behavior of the shear viscosity ri(T) [|24| , P6t |28|| . OTP is a 
relatively simple compact, non-polar molecule with an overall shape that comes close 
to the spherical particles preferred by theory. It consists of small, hopefully rigid 
molecular units which interact via weak non-directional van der Waals forces. Therefore, 
we expect that the dynamics of supercooled OTP represents structure independent, 
generic properties of glass forming liquids and does not depend on microscopic structural 
peculiarities. 

The molecule consists of a central benzene ring and two lateral phenyl rings in 
ortho position. For steric reasons, the lateral two phenyl rings are necessarily rotated 
out of plane as shown in figure |l]. In addition, the overcrowding in the molecule leads 
to significant bond-angle and out-of-plane distortions of the phenyl-phenyl bonds. Such 
structural irregularities may explain why OTP can be supercooled far easier than m- or 



p-terphenyl p9|. 



A number of structural features of the OTP molecule have their origin in the close 
proximity of the two bulky lateral phenyl rings. However, the details of the molecular 
structure may depend on whether it is in the crystalline, liquid, glassy or gas phase. 

In the crystal, the angles for the out-of-plane rotation of the lateral phenyl rings 
are 0i ~ 43° and 02 — 62° |^0|, For isolated molecules, an old electron diffraction 



study had suggested 4>i = (p2 = 90° but newer experiments and calculations agree 




Figure 1. Molecular geometry of orthoterphenyl (C18H14) as observed in the crystal 
structure. It is constituted by three phenyl rings, the two lateral rings being attached 
to the central one by covalent bonds. 



that in the gas or hquid phase 40° 



< 



< 65° [33, 



35|. 



The crystal structure belongs to the orthorhombic space group P2i2i2i with four 
molecules per unit cell and lattice parameters a = 18.583 A, b = 6.024 A, c = 11.728 A 

A sketch of the structure is given in figure 1 of 



a 31 



at room temperature 
Ref. Experimental studies of the crystalline state have been performed mainly 

on polycrystalline powder samples or on rather small single crystal needles. However, 
single crystals of high quality and considerable size (several cm^) can be grown out of a 
hot methanol solution either by very slow cooling or by evaporation over several months 
36| . The crystals grow preferentially along the shortest axis b. 

Crystallization behavior of high purity OTP has been studied extensively over a 
range of supercooling up to above the melting point in connection with the glass forming 
ability fS^, |8|, ||, 



41 



For many glass formers the rate of cooling is taken as 
determining whether or not a liquid will form a glass. If unperturbed, OTP and some 
liquids of similarly shaped molecules require years to nucleate |37, However, if 
abundant stable nuclei are artificially created in the sample, then the overall liquid 
solid transformation occurs fairly rapidly on the time scale of seconds or minutes. The 
crystal growth rate of OTP at atmospheric pressure passes through a maximum value 
at 312 K |3^. On the other hand, nucleation frequencies are difficult to measure and 
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are usually calculated using classical nucleation theory. For example, the steady state 
homogeneous nucleation frequency is highest at 251 K The bulk heterogeneous 

nucleation frequency maximum lies between 260 K and 275 K for contact angles of 100° 
and 70°, respectively [0. In most cases, nucleation takes place (heterogeneously) on 
surface heterogeneities whereas homogeneous nucleation is rather exceptional in glass 
forming liquids [0. Heterogeneous nucleation at glass surfaces has not been studied 
rigorously, even though it is the most common origin of devitrification. Such nucleation 
may occur whether the surface of the glass is in contact with a container wall or exposed 
to atmosphere. Heterogeneities characterized by small contact angles have a marked 
effect on the cooling rates required for glass formation, with the required cooling rate 



increasing with decreasing contact angle Nevertheless, with some care it is 



relatively easy to supercool OTP below the melting point. Depending on the cooling 
rate, the caloric glass transition temperature Tg lies between 242 K and 245 K |2^ . 
Thus, OTP is rather easy to handle as the melting and the glass transition temperature 
are close to ambient temperature. 

One further advantage particularly useful in neutron scattering studies is that the 
molecule can be selelctively of even fully deuterated. As a proton containing system 
the scattering is almost purely incoherent enabling us to study the single particle or 
tagged-particle dynamics. Fully deuterated the scattering is mainly coherent giving rise 
to the observation of the pair correlation function. Partial deuteration can be used to 
hide groups of atoms and to study the intra- and intermolecular motions separately. 

Furthermore, we chose OTP for our investigations since it has been extensively 
studied in the literature by many spectroscopic methods as we will see in the next 
section (|^). Thus, as a function of temperature and pressure many data are available 
for comparison. 
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3. Glass Transition in OTP 

3.1. Macroscopic Quantities 

The calorimetric techniques used for measuring the heat capacity cp of glasses are the 
same as used for other solid materials. In particular, differential scanning calorimetry 
(DSC) |]43| , or differential thermal analysis (DTA) |^ are used to determine the 
glass transition temperature Tg. When the heat capacity cp(T) of a supercooled liquid 
is measured as a function of temperature (or pressure) one notices a significant change 
occurring at a temperature Tg (figure |^). A similar effect is observed in volume- 




ICO 
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Figure 2. Specific heat of liquid, supercooled liquid, glassy and crystalline OTP 
obtained by differential scanning calorimetry. The data in the left part (a) are obtained 
with a scanning speed of 0.3 K/s (reproduced from |Q) while for the data on the right 
(b) the scanning speed was 0.003K/s (values taken from |4^). 



temperature-pressure measurements reproduced in figure ^. A "discontinuity" in 
the thermal expansion coefficient a = dlnV/dT\p, the derivative of the volume V at 
constant pressure P can be deduced from figure ^ for the glass while nothing particular 
happens in the crystal. Although these features are characteristic of a second-order 
equilibrium phase transition, Tg cannot be related unambiguously to an equilibrium 
function, such as the free energy. In particular, the position and the sharpness of 
the transition near Tg depends on the cooling rate at which for example the heat 



are 



capacity is monitored. In figure ^ results with two different cooling rates ^ 
compared. Conventionally, cooling rates of about lOK/min are used to determine Tg. 
Operationally, Tg is defined as the center of the temperature range of the large change 
in the heat capacity in a DSC run with cooling rate T. Note, that in this definition 
the glass transition is related with a low-frequency excitation: e. g. for a typical rate 
of T = lOK/min ~ 0.2K/s and a width of the cp-jump AT = lOK the excitation 
frequency is = T/AT = 0.02 s~^, thus probing correlation times of about 50 s. The 
glass transition temperature defined in this manner can be understood as an isodynamic 
point which allows for a comparison of different glasses. 




Figure 3. Varations of the isobaric specific volume of liquid, supercooled, glassy 
and crystalline OTP as a function of temperature for various pressures 0.1 MPa (□), 
50MPa (■) and 100 MPa ({>)■ The solid lines are linear extrapolations from either 
side. Reproduced using polynomial coefficients given in ||4^ . 



When pressure is applied the transition point shifts to higher temperatures as can 
be seen in the volume-temperature-pressure relation in figure ^. Note, that under 
pressure OTP shows a very strong density change. Calorimetrically, the glass transition 
temperature in OTP as a function of pressure, Tg(P), has been measured by DTA 
by Atake and Angell Up to at least about 100 MPa the variation of Tg{P) is 

linear with a slope dTg/dP ~ 0.26K/MPa (260K/GPa). These values for OTP and 
its isomeric mixtures are the highest yet observed for molecular liquids (though some 
polymers exhibit higher values). It was suspected that dTg/dP is larger the more closely 
the molecular substance behaves repulsively like a hard sphere. These facts have been 
further motivations for us to use OTP as a model system also for pressure dependent 
neutron scattering experiments (section 



3.2. Primary a-Relaxation: Stretching and Scaling 

The fiuctuation-dissipation theorem connects a linear response with the time evolution 
or spectral distribution of spontaneous fiuctuations. Thus, relaxation in the equilibrium 
state can be studied by a variety of spectroscopic or scattering methods. 

Time-resolved measurements are performed for instance by photon-correlation 
spectroscopy (PCS) which measures ultimately the evolution of density correlations 
$(t) ^ The PCS data in figure || illustrate nicely the most prominent features 



of structural relaxation: stretching, scaling and a strong temperature and pressure 
dependence of the mean relaxation time. 
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Figure 4. Normalized time-correlation function of liquid OTP, obtained by photon 
correlation spectroscopy in depolarized (VH) scattering geometry at constant pressure 
(top) and at constant temperature (bottom). The Kohlrausch stretched exponential 
is fitted to the data (dots) with exponents (3 between 0.53 and 0.62. Reproduced from 
For comparison, a simple exponential function was added (long dashed line). 



The decay of the time correlations $(t) is called stretched since it extends over a 
much wider time range than a simple exponential process. An exponential decay is 
expected when the dynamics consists of uncorrelated random steps. This is the case 



for tagged-particle diffusion like motion in the long time and long distance limit |50 



On microscopic scales, however, the motion of molecules is a highly collective process 
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marked by memory effects. 

In many cases the structural relaxation function is well described by a 

Kohlrausch law |^T], obtained by replacing t/r in the Debye law by (t/r)^: 

<l>(t)=exp(-(t/rf), 0<P<1. (1) 

For this function, the mean relaxation time (r) can be calculated as 

POO 

{r) = / dtm = rr'T{(3-') (2) 
Jo 

where r() is the gamma function. The exponent /3 is not (!) a material constant; it 
differs for different correlation functions or observables which are probed by different 
experimental techniques [Q (see also figure ^. Such variations of /3 demonstrate that 
an understanding of structural relaxation will not be possible without accounting for 
the different ways in which experimental probes couple to the microscopic degrees of 
freedom. 

On cooling or compressing a viscous liquid, (r) increases rapidly as can be seen in 
figures § and ^ while the line shape of $(t) , expressed through f3 varies only weakly. This 
scaling property, known as the time-temperature-pressure superposition, implies that in 
a first approximation the measured individual correlations $(t; T, P) can be collapsed 
onto a master curve 

<!>{t-T,P) = <l>{t/T{T,P)) (3) 

by rescaling the time axis to i = t/riT, P). This property turns out to be particularly 
useful in the analysis of neutron scattering experiments, to be discussed later in 



sections 3.5 and 7.5. 



In many spectroscopic experiments, the response of a viscous liquid is probed as a 
function of frequency. A complex susceptibility x(^) is related to the relaxation function 
$(t) through 

/■oo 

XH = Xoo + iXoo - Xo) / dte'-'m. (4) 

Jo 

For high frequencies, the system reacts as a solid x(^ ^ t^^) = Xoo] for very low 
frequencies the fluid-like response is given by x(i^ ^ t^^) = Xo- 

For intermediate frequencies u ~ the Kramers-Kronig relation relates the step 
in the real part x'{^) to the maximum in the dissipative part x"i^)- An exponential 
decay of $(t) leads to 

= Xoo + . (5) 

1 — lUJT 

which in the context of structural relaxation is often named after Debye, who derived 
it for Brownian motion of uncorrelated dipoles in a viscous medium . 



For a stretched relaxation function, the step in x'{^) ciiid the peak in x"{^) ^ire 
much broader than a Debye peak. To describe x{^) one can either decorate (^ with 
one or two exponents ^ ^ or use the Fourier transform of (|I|) 
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Frequency dependent dielectric susceptibilities for pure OTP have rarely been 
studied ^ 0, |6^ probably due to the extremely weak dipole moment. 
Instead, in many studies the reorientation of polar host molecules is measured 
p^ , |65| , [66| , |67| , |6^ , |69| , [70| , [7I[] . Typical temperature dependent dielectric spectra and 
the corresponding master curve of 4.2 weight-% anthrone in OTP are shown in figure ^. 
For comparison, a Debye peak is also shown. Molecular reorientation of different probe 
molecules has also been studied by optical |73l, ESR 0, |75| and NMR spectroscopy 
7|. 




Figure 5. (a) Dielectric loss of 4.2 weiglit-% anthrone in OTP. Data are scanned 
from [|64[ and re-fitted with a Fourier transformed Kohlrausch function. The exponent 
P seems to increase systematically from 0.52 at 255.5 K to 0.57 at 267.1 K. — (b) 
Nevertheless, in a good approximation, all data can be condensed onto a master curve. 
Same symbols as in (a). For comparison, the dashed line shows a Debye peak. 



In figure ^ some published values of the Kohlrausch stretching exponent f3 are 
compared. For a given experimental probe, (3 increases systematically with temperature 
in this temperature range. Thus, the time-temperature superposition (j^) only holds in 
a first approximation and in a limited temperature range. In addition, for any given 
temperature, different experiments see different exponents (3 although approximate 
values of (3 seem to be correlated with other characteristics of the glass transition 
dynamics JTHI, |79 . 



3.3. Temperature and Pressure Dependence of the a-Relaxation Time 

We have seen that the characteristic time (or frequency) of structural relaxation strongly 
depends on temperature and pressure. Cooling OTP by only 5K can slow down the 
decay of the correlation function by a factor of ten or more, as shown in figures ^ and |. 

Figure ^ demonstrates that over many decades the characteristic time r observed 
by different spectroscopic techniques follows the same temperature and pressure 
dependence, although at a given temperature or pressure the absolute values of r may 
differ. This finding is confirmed by measurements on many other glass formers. 
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Figure 6. Stretching exponent (3 from Kohlrausch fits to tliermal effusion (x, OTP 
witli 9% o-phenyl-phenol j?^), plioton correlation (♦ |47j), and dielectric loss (□, 
4.2% anthrone in OTP, our fits to the data ^ of Fig. |). 



The famous Maxwell relation rj = G'oo(t) connects the shear viscosity r] with the 
relaxation time over the shear modulus at infinite frequency Goo- Under the assumption 
that Goo varies only weakly with external control parameters like temperature and 
pressure, it follows rj oc (r). 

Therefore, another operational definition of the glass transition temperature Tg is 
given via the temperature, where the viscosity rj reaches ^^(Tg) = lO^^Pa-s. With shear 
moduli of solids, typically of the order Goo — 10 GPa, this corresponds to a relaxation 
time of the order of 100 s. This is the same time scale on that structural relaxation falls 
out of equilibrium in a caloric experiment (see above). 

The Maxwell relation makes it possible to determine the temperature or pressure 
dependence of the structural relaxation time beyond the frequency range of any 
spectroscopic technique. In fact, the viscosity of OTP has been measured over more 



than 14 decades as a function of temperature and pressure |]2^, |25|, The 



temperature dependent data at ambient pressure are shown in figure |^ together with 
viscosities of other glass formers in a log// vs. Tg/T representation (Angell plot) 
motivated by the naive expectation of thermally activated transport ||8^. Although 
never observed in any real glass former the Arrhenius law r] oc exp{—EA/T) is taken as 
a reference to classify the temperature dependence of ri(T, P). 

This classification according to Angell is physical since the curvature in the 
Arrhenius plot is correlated with the structure of a glass former: strong liquids, which 
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Figure 7. Temperature dependence of the isobaric relaxation time of viscous OTP for 
pressures between 0.1 MPa and 125 MPa from PCS (open symbols from right to left in 
steps of 25 MPa) ]48| a nd dielectric spectroscopy (closed symbols from right to left in 
steps of 19.6MPa)|§. 



deviate only slightly from the Arrhenius law are typically networks like Si02 or Ge02 
where the transport is controlled by thermally activated breaking of covalent bonds. 
Largest deviations are found for molecular liquids, ionic melts and some polymers. They 
are termed fragile liquids. In this classification OTP is one of the most fragile liquids, 
which is another motivation for its recurrent use as a model system. 

From free-volume or entropy arguments it is possible to derive the popular 
Vogel-Fulcher-Tamman fit formula 

7]{T) oc exp ( ^ ) (6) 



T-Tr 



with some arbitrary reference temperature Tq <^ Tg. In any case, (|^) and similar 
functions are not able to fit f]{T) data over the entire temperature range, but they 
work much better when applied to a restricted range. A critical review of these and 
other approaches can be found in Ref. |^. The differential analysis of ri{T) in the 
inset of figure ^is one way to show that the temperature dependence of ri(T) cannot be 
described by a single physical model over the whole temperature range. In particular, 
one recognizes a change of regime at about 290 K. Within mode-coupling theory this 
is interpreted as a cross-over from liquid-like to hopping dynamics (section For 
OTP the value ^ 290 K was first determined by strict elastic neutron scattering 
18^ (section |6.3|) and in the meantime it reappeared in many other measurements and 
simulations, as will be shown throughout this review. In its simplest version, mode- 
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Figure 8. "Angell plot" : Viscosities of various glass forming liquids in an Arrhenius 
representation with reduced temperature scale Tg/T. This plot motivates Angells' 
"strong" — "fragile" classification in which OTP is one of the most fragile liquids. 
Reproduced from ]2^ . 



coupling theory predicts for T > a, power law divergence 

r/ oc (T - T^y, 1.76 < 7 < 4.83. (7) 

Figure D shows a three parameter fit of r]{T)/T yielding a critical temperature = 290 K 
and an exponent 7 = 2.59. 

3.4- Diffusion and Decoupling of Time Scales 

It is interesting to compare the translational and rotational diffusion coefficients, Dt and 
Dr, which are integrals over single particle correlation functions, with the shear viscosity 
rj. The shear viscosity depends upon intermolecular momentum transfer through cross 
correlation functions of neighbouring molecules. In simple van der Waals liquids, the 
validity of the Stokes-Einstein and the Debye relations, 

A=^, Dr = ^, (8) 

b-nrjrs aTrrjri) 




Figure 9. Temperature variation of ri{T)/T. The solid line is a free fit of a power 
law (T — Tc)^ expected from mode coupling theory yielding Tc ~ 290 K and 7 = 2.59. 
The arrows indicate the caloric glass transition temperature Tg — 243 K, the critical 
temperature Tc = 290 K and the melting temperature = 329 K. In the inset 
(taken from ]6^ ) the temperature dependence of the maximum fmax of the dielectric 
a-process (o) and the viscosity ri{T) (x) is plotted as [— dlog(a;)/d(l/T)]. Such a 
representation linearizes a VFT-type temperature dependence and an Arrhenius 
behavior appears as a horizontal line. Tb ~ 290 K marks the temperature where the 
data begin to deviate from a VFT-fit (line). Above Ta — 455 K ri{T) follows an 
Arrhenius law. 



originally derived for a Brownian particle in a viscous liquid, is well established in 
the fluid regime, although it is a non-trivial fact. The lengths rs = 0.21 nm 
~ r/) = 0.23 nm ^ are temperature independent and smaller than the geometrical 

0.37nm 



van der Waals radius rw = 0.37 nm It is most remarkable that Dt and Dr are 

proportional to ?7~^ in the whole fluid regime. 

Rotational self diffusion in neat OTP has been measured by ^H-NMR stimulated 
echos translational self diffusion by ^H-NMR static gradient stimulated echos 

3^, |88|. Towards lower temperatures, tracer diffusion coefficients of two photochromic 



tracer molecules having a size comparable to OTP and dissolved in OTP at low 
concentrations were measured by forced Rayleigh scattering |8^. Results are shown 



in figure 10 



Rotational self diffusion could be measured almost down to Tg where Debyes' 
relation (cf. still holds. Apparently, momentum transport is closely connected with 
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Figure 10. Log-log plot of the translational and rotational diffusion coefficients 
Dt and Dr of OTP vs. viscosity r/ and temperature T: translational self diffusion 
(•); translational tracer diffusion (different photochromic molecules in OTP, A,n); 
rotational diffusion coefficient (♦,0) scaled along Dt/Dj- — 4/3r^ (see (^) using r = 
0.22 nm). The steep line represents the Stokes-Einstein relation Dt = k-^T / [Q-K-qr s) 
using rs — 0.22 nm and fits well the I?t(^) data above Tc ~ 290 K and the Dr{f]) data. 
The flat line that fits the Dt{r]) below Tc ~ 290 K has a slope 7^~o '^^ and crosses the 
77^^ line at 290 K. Reproduced from |]90| . 



molecular reorient at ions, r Translational diffusion, however, decouples from rotational 
diffusion and viscosity and crosses over around = 290 ± 10 K to a weaker temperature 
dependence approximately given by Dt oc 7/^° ''^. Thus, Dt and r] behave differently as 
the mechanism of molecular motion changes from fluid to glass like behaviour. On the 
other hand, D^ remains proportional to rj"^ at temperatures down to Tg indicating 
that rotational motion remains more efficiently coupled to structural relaxation on 
approaching the glass transition, at least, in fragile glass forming liquids like OTP, 
toluene ||9T|, salol, glycerin and others |92 . 

Indications for a cross-over between two distinct dynamic regimes at a temperature 
Tc above the glass transition temperature Tg have been found for many simple 
supercooled liquids |9^. A discontinuity for the temperature coefficient of r/?7 or {Dri)~^ 
demonstrates the failure of the Stokes-Einstein relations for T < T^. The existence of 
the cross-over temperature Tc is further substantiated by a plot of the correlation times 
on a reduced temperature scale Tg/T, which then shows a bifurcation of the two relevant 
a and /? processes for T < Tc. 
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4. Experimental Method 
4.I. Neutron Scattering 

In the previous section we have presented a variety of experimental glass transition 
signatures observed with conventional spectroscopic methods. What can we expect to 
learn about the glass transition from neutron scattering? 

The energy of neutrons of the appropriate wavelength for structural and dynamic 
studies corresponds to thermal energies for temperatures between a few degrees K to 
well above room temperature. In fact, neutrons produced from a reactor or spallation 
source have their wavelength "adjusted" by the equilibration at a given temperature; 
thus "cold" neutrons have a wavelength of several A, "hot" neutrons have a wavelength 
of a fraction of an A and "thermal" neutrons have a wavelength close to 1 A. Because of 
this wavelength-energy link, neutron inelastic scattering allows the simultaneous study 
of spatial structure, on the one hand, and dynamics, on the other, at the atomic and 
molecular level 0, ^, 0, 0. 



Neutron scattering experiments on viscous liquids are motivated by the following 
observations: In amorphous and crystalline solids vibrational excitation frequencies are 
generally of the order of 1 THz; long wavelength phonon frequencies of sound waves are 



found in the GHz to THz region depending on the wave number Q [^, |9^. Furthermore, 
the mean relaxation time in normal fluids is typically of the order of Ips ||9^. Thus, 
for the observation of the onset of the glass transition - when relaxational motions 
decouple from vibrational ones - we have to look at the GHz to THz range which is 
traditionally the domain of scattering experiments. Such high frequencies are nowadays 
accessible with many experimental techniques like e. g. dynamic light scattering, inelastic 
X-ray scattering, quasielastic neutron scattering, dielectric spectroscopy or impulsive 
stimulated light scattering (section ^731) . However, not all of these techniques have been 
applied to OTP. 

Neutron scattering remains one of the most versatile methods for the study of the 
fast dynamics. It can be used in almost every type of material since it does not depend on 



optical transparency, dipole moments or polarizabilities |9^. In neutron scattering one 



observes excitations at a well defined momentum transfer Q of the order of several A~^, 
i. e. microscopic motion on interatomic and intermolecular distances can be resolved. 

The observed quantity is the so called dynamic structure factor S{Q, lj) and is a 
well defined correlation function [0, 

^(Q,.;) = 1 /" ^e^-* (e.QR,(o)g-iQR,(t)^^ (9) 

^ j,k=l 

which couples directly to the nuclear coordinates Rj(t). It therfore provides an ideal 
link between experiment and theory as S'(Q, uj) is the space-time Fourier transformation 
of the well know van Hove correlation function G{r,t) [^]. 

Equation (H) is based on the coherent superposition of partial waves scattered by 
pairs of nuclei j, k. The coherence can be diminished by spin flips except when both 
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Table 1. Neutron scattering cross sections for incoherent and coherent scattering and 
absorption cross section in barn ''— 1on^•T>^2^ 



atom / molecule 




fcoh 


CTabs 


H 


80.27 


1.76 


0.33 


D 


2.05 


5.59 


0.00 


C 


0.00 


5.55 


0.00 


OTP-rfo 


1123.8 


124.5 


4.62 


OTP-dio 


341.6 


162.8 


1.32 


0TP-di4 


28.7 


178.2 


0.00 



(12) 



partial waves are scattered by the same nucleus j = k. This leads to an additional 
incoherent scattering 

which contains valuable information about tagged-particle motion. 

In an actual experiment one measures always the sum of both contributions 

S'exp(Q, Uj) = Crcoh5'coh(Q, t^) + 0'mc5'inc(Q, t^)- (H) 

Incoherent and coherent scattering can be separated by a polarization analysis |9^, |100|| 
which is based on the different scattering cross sections for spin-flip a-^ = 2/3ainc and 
non-spin-fiip cr|| = Ccoh + ^/^(^inc 

S^^iQ,^) = Crcoh'S'coh(Q7 ^) + |crinc5'inc(Q5 '^)5 
'S'ti(Q)^) = 30"inc'S'inc(Q, ^^)- 

Although this implies a considerable loss in neutron flux in high resolution experiments 
such experiments are indeed feasible as shown below. 

In practice, most often one concentrates on either incoherent or coherent scattering 
by choosing appropriate samples. In organic samples like OTP, the total cross section 
is largely dominated by incoherent scattering from hydrogen (table [^). Mainly coherent 
scattering can be measured on perdeuterated OTP-d^. Partial deuteration has been 
employed to hide the lateral phenyl rings |p.01|| (section |0|) . 

The spectrometers used for neutron scattering on OTP are listed in table 
They belong to basically five types ||103|| and are all located at the Institute Laue- 
Langevin in Grenoble, France. Three of them measure S'(Q,ti;), one the so called 
intermediate scattering function S{Q,t) which is the Fourier transform of S'(Q,u;) and 
one spectrometer works without energy analysis: 

• Triple-axis spectrometers (IN12) are instruments which can measure the four 
dimensional dynamic structure factor ^(Q,^^). In single crystals the phonon 
dispersions relations uj{Q) can explicitly be measured. 



Neutron Scattering Studies of the Model Glass Former Ortho-terphenyl 



19 



• Time-of-flight (TOF) spectrometers with multichopper (INS) or crystal (IN6, D7) 
monochromators where the energy of the scattered neutrons is measured by the time 
they need to travel a certain distance are used to measure vibrational dynamics as 
well as fast relaxation. 

• Back-scattering (BS) spectrometers achieve high energy resolution by (nearly) 180° 
Bragg scattering from both the monochromator and the analyzer crystals. They 
are used to measure the amplitude of elastic scattering and the onset of the slow 
relaxation. In order to measure full spectra the incident energy is changed either 
by Doppler shifting the neutron wavelength (INIO, IN16) or by thermal expansion 
of the monochromators (INIO, IN13). 

• Neutron spin-echo (NSE) instruments (INll) require by construction the use of 
polarized neutrons. They profit from the fact that the magnetic moment of the 
neutron undergoes Larmor precessions in a homogeneous magnetic field. This is 
used to encode the neutrons' energy in the number of Larmor precessions. The 
polarization analysis of the scattered neutrons amounts to an implicit Fourier 
transform of the scattering function S{Q, u) and has a much larger dynamic range 
as frequency domain spectrometers at the cost of a lower contrast. 

• Diffractometers (D7, D20) are used to measure the static structure factor of hquids 
and amorphous solids. The integration over all energies is intrinsically performed. 
In addition to a time-of-flight option, D7 can perform a polarization analysis and 
is thus able to measure simultaneously incoherent and coherent scattering. 

For isotropic samples like powders, liquids and glasses we can only access Q = |Q|. 



4-2. Data Evaluation and Limitations 

To obtain the desired quantity S{Q, u) from the measured scattering intensities several 
standard corrections have to be performed like normalization to the incoming flux 
and to the scattering of a vanadium standard, subtraction of container scattering and 
background noise and, if necessary, calculation of attenuation and absorption effects and 
conversion to energy. 

Since all instruments measure at given scattering angles 6 a conversion to the 



physically interpretable wave number Q has to be done. In figure [Tl| we show the 
dependence of the wave number Q on the energy transfer u for several scattering angles 
29 = . . . 180°. For backscattering and neutron spin-echo instruments, where the energy 
Ei of the incident neutrons is much larger than the energy transfer u the momentum 
transfer Q at a given scattering angle is essentially constant (left panel) . However, when 
the energy transfer uj is of the order of Ei - for instance in time-of-fiight spectroscopy 
- then the dependence of Q on cannot be neglected (right panel) although sometimes 
still done PI . 



In any substance, one observes at intermediate temperatures significant 
contributions to S{26,uj) up to several lOmeV. Before calculating a Fourier transform 
to get S{Q,t), S{29,uj) has to be interpolated to S{Q,lj). As illustrated in figure 0, 
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Table 2. Neutron spectometers used in our studies of the statics and dynamics of 
OTP summarized in this review. All instruments arc located at the Institute Laue- 
Langevin in Grenoble and use cold neutrons except IN13 which uses thermal neutrons. 
Only IN12 is a single detector instument, on all other instruments use multidetector 
banks which cover a wide range of scattering angles. The energy resolution 5E varies 
with detector angle and chosen incoming wavelength Aj. Often quasiclastic scattering 
data are analyzed by Fourier deconvolution and the approximative time ranges are 
given. The spin-echo spetrometer INll measures directly in the time regime and can 
be operated with a single- or a multidetector. 



Instrument 


type 


Xi (A) 


5E (/xeV) 


Fourier time (psec) 


INS 


time-of-flight 


5.7-7.5 


25-40 


0.2 - 35 


IN6 


time-of-flight 


5.1 


80-120 


0.1 10 


INIO 


back-scattering 


6.27 


1-3 


200 - 2000 


INll 


spin-echo 


5.5-6.0 




1 - 2200 


INI 2 


triple-axis 


1.55 (kf) 


300 




IN13 


back -scattering 


2.23 


8-10 


10 - 200 


IN16 


back-scattering 


6.27 


1 


100 - 2000 


D7 


time-of-flight 


4.8 


150 


0.2-5 


D7 


diffractometer 


4.8 






D20 


diffractometer 


0.82 








Figure 11. Maximum dynamic range of neutron scattering: (a) on a logarithmic 
energy scale to emphasis the small energy transfer region typical for BS or NSE 
instruments and (b) on a linear energy scale to emphasis the large energy transfers 
typical for TOF instruments. The lines display the frequency dependence of the wave 
number Q{26,w) for several scattering angles ranging form 0° to 180° in steps of 20° 
from bottom to top. In any real experiment the subregion depends on the largest and 
lowest scattering angle accessible by each individual instrument. Only the data along 
aQ = const-line in the gray hatched subregion can be used when interpolated spectra 
S{Q, to) arc Fourier transformed. The vertical line a; = is the lower energy bound 
as the elastic channel is needed for the FT. At a given wave number Q the highest 
accessible energy is restricted and limited by the lowest scattering angle 9min- 
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this imposes restrictions upon the accessible energy transfers. The Fourier transform 
S{Q, t) can only be evaluated if the frequency spectrum includes the elastic channel 
= 0. The truncation problem at ci^max is mathematically solved by a multiplication 
with the step function and leads to a lower bound in Fourier time t below that the 
Fourier transform is not reliable ||105|| . The truncation effect may also appear in 



S{Q,t = 0) = 2 duS{Q,u) which is then less than the true static structure factor 
S{Q). Calculating the density and tagged-particle correlation function ^{Q, t) from the 
dynamic structure factors 

$(g,t) = 5,oh(g,t)/5coh(Q) and <l>'{Q,t) = S,UQ,t)- (13) 

might then lead to a wrong normalization. For incoherent scattering this problem can be 
circumvented by measuring S{Q,uj) at low temperatures where the scattering is totally 
elastic and S{Q) = 1. For coherent scattering this is not possible as S{Q) strongly 
depends on temperature and pressure. Note, that the truncation effect depends on Q 
and on the dynamics of the studied system itself and that two competing effects are 
operative: With decreasing Q the accessible energy range gets narrower (see figure |TT|), 
but at the same time the quasielastic broadening and the inelastic intensity decrease 
as well. Though, while it is straightforward to measure a plateau or a decay of the 
correlation function it is rather difficult to determine its absolute value within a few 
percent. The truncation effect can be studied using several upper cut-off frequencies 
and is in general negligible at least for intermediate and large wave numbers Q. 

The determination of the dynamic structure factor in absolute units is further 
complicated by the inevitable presence of multiple scattering. Multiple scattering 
is basically a convolution of the dynamic structure factor S{Q, u) with itself. The 
correction for it is far from being a standard data treatment. In the inelastic spectrum, 
most multiple scattering intensity comes from elastic-inelastic scattering histories [[L06|] . 
This contribution is nearly isotropic and therefore weakest at large scattering angles. It 
tends to smear out the characteristic Q dependence of the scattering, and dominates at 
small angles where the single-scattering signal is expected to vanish as Q^. 

There are several different approaches to correct for this effect. One can perform 
a Monte-Carlo simulation of the neutron histories in an effective way [|107|| . In order 
to calculate the fraction of multiply scattered neutrons at one point {Q, u) one has, in 
principle, to know the ideal scattering in the whole {Q, c<j)-space in absolute units. As 
can be seen in figure |lT] only part of it is experimentally covered. Therefore, on needs 
a physical model for the ideal scattering including Q- and temperature dependences. 
This makes the correction even more complicated for coherent scattering where elastic, 
quasielastic and inelastic scattering are oscillating functions of Q. 

A new experimental, model independent, however very time consuming approach 
for multiple scattering correction in inelastic neutron scattering consists in the use of a 
set of different incident neutron wavelengths in connection with Monte Carlo simulations 
108||. Combining data at different wavelengths one can reliably establish an ideal, single 



scattering function in a wide Q range which then serves as an input to the simulation 
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followed by a self-consistent, iterative correction. 

Another possibility at least for incoherent scattering is to expand the scattering 
function S{Q,u!) = A{u) + Q'^Bioj) + . . . which enables us to estimate the multiple 
scattering Aioj) by — extrapolation from the S{26{Q,u}),uj) spectra however, 
loosing the elastic channels |36[. In practice, one often uses thin samples with high 
transmissions of the order of 90% and more. Then multiple scattering is in general 
neglected at least at high scattering angles [|107| . 

Recently, an alternative approach has been worked out. By extensive Monte Carlo 
simulations on very simple models some generic trends of multiple scattering were 
identified [|109|| : it distorts the wave number dependence much more than the frequency 
dependence of S{Q,uj). 

Neutron scattering is always a compromise between the desired statistics, the wave 
number resolution and coverage, the elastic energy resolution and the accessible energy 
window. Apart from notoriously lacking statistics, one major problem in neutron 
scattering is that practically all high resolution spectrometers lack dynamic range. It 
is known that in order to measure the dynamics of supercooled liquids and glasses 
properly one has to cover several decades in frequency or time. Several approaches 
(partly) overcome this problem: (i) the use of master functions obtained by scaling 
procedures like the time-temperature superposition and / or (ii) the combination of 
data from different spectrometers with overlapping time ranges. Over the years we have 
been using both approaches. Elastic resolution effects can be corrected in the S{Q, t) 
representation by division with the Fourier transformed resolution function (Fourier 
deconvolution) that has been measured at very low temperatures either on the sample 
itself or on a Vanadium standard. 

Neutron spin-echo spectrometers which measure directly in time domain have a 
rather large dynamic range of about 3 decades and more, depending on the scattering 
angle and the incoming wavelength [[103|] . However, the concomitant shortcuts are low 
countrates and low Q resolution compared to multidetector instruments. Recently, a 
multidetector (42 detectors) has been installed on the NSE instrument INll covering 
about 40°. Now, by binning several detectors at least the statistical quahty can be 
improved. The Q resolution is still implicitly limited by the transmission band of the 
velocity selector of typically AA/A ~ 10 % leading directly to AQ/Q > 10%. One 
should be aware that in such a case one measures averages ($(Q,t))Q over correlation 
functions. Increasing the Q resolution by a narrower band width using for example a 
graphite monochromator leads to an unacceptable reduction in intensity. 

Improvements of neutron sources [|110| , |111|| to increase the incoming flux and 
of instrumental performance ||112|| to increase the efficiency are currently under way. 
In order to test the new developments and extentions of the theory (section |5.4| ) 
experimentally or to attack even new questions the dynamic range and concomitantly 
the Q-resolution and the Q-coverage have to be enhanced. 

Let us demonstrate to which extent it is desirable to push improvements for a 
real progress. As an example we show in figure |l^ the tagged-particle and the density 
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dynamic structure factor of OTP measured in the only proper way by polarization 
analysis on the instrument D7 at the ILL. Currently, D7 is the only multidetcctor 
instrument in the world which is able to separate the energy resolved incoherent and 
coherent scattering from one sample. Some triple-axis spectrometers can do polarization 
analysis as well, however, are equipped with a single detector only. On each other 
instrument one always measures both contributions which are differently weighted 
depending on the spectrometer type. This possibility is particularly interesting for some 




1 
CO (meV) 



1 2 

CO (meV) 



Figure 12. (*) Coherent (open symbols) and incoherent (full symbols) dynamical 
structure factor S{Q,u}) of OTP for several values of Q at 340 K as measured on D7 
with polarization analysis. Note the different quasielastic scattering for different Q. 

important questions not only in glass science which have not been studied up to now. For 
example, it would be possible to study exclusively the center-of-mass (COM) motion 
of OTP. Measuring the incoherent scattering of OTP-dio where only the lateral phenyl 
rings are deuterated gives access to the tagged-particle correlation function of the COM 
motion. Conversely, measuring the coherent scattering of central ring deuterated OTP- 
(^4 would lead to the density correlation function of the COM motion. At the same time, 
important information can be gained from the coherent scattering and the incoherent 
scattering of the lateral phenyl rings in OTP-rfio and 0TP-(i4, respectively, about e. g. 
relaxational modes and their time scales or vibrational excitations of intramolecular 
motions. 

Furthermore, it would also be possible to compare both structural relaxation times 
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and their Q dependence in the context of reccurent used models hke "de Gennes" 
narrowing ||113|| . The advantage is evident: such studies could be performed on one 
and the same sample under exact identical experimental conditions. 



The spectra shown in figure |T2| were obtained under fairly optimized experimental 
conditions. Nevertheless, although the transmission of 67% was extremely low and 
the energy resolution rather moderate data accumulation took about 42 h! Such test 
experiments make it clear that, for example, increasing the incoming flux by a factor 
of two while leaving the instruments as they are is far too less. The total efficiency 
which is a product of the source strength, the transport efficiency and the instrument 
efficiency has rather to be increased by orders of magnitudes [|112|| . 



4.3. Other Experimental Approaches 

In section ^ we have already presented many experimental techniques probing the glass 
transition phenomena and the slow structural relaxation of OTP. Therefore, we limit 
ourselfs to a short enummeration of those other experimental approaches that probe the 
fast dynamics in the GHz - THz band and have been applied to OTP. 

Light scattering spectroscopy has been often employed in OTP and many other 
supercooled liquids and glasses. Raman scattering has primarily been performed in 
connection with the Boson peak and the study of the vibrational dynamics of the glassy 
and polycrystalline state ||1 14| , |115| , |116|| . Brillouin scattering explores the longitudinal 
(LA) and transverse (TA) acoustic modes in the GHz range and their interaction with 
the structural relaxation [[IIT] , |118| , |119|| . Inelastic X-ray scattering extends the Q-lo 
range of the conventional Brillouin technique and makes it possible to study sound 
waves on intermolecular distances and in the THz range ||120| , |121| , |122|| . Depolarised 
light scattering using geometries in which the usual LA and TA Brillouin components 
are forbidden is the most often used method to study the sub GHz to THz dynamics 
|123| , |124| , |125| , |126| , |127|| . Although the scattering mechanism and the connection 



to density fluctuations are still debated ||128|| the experimental results have often been 
compared to mode-coupling theory. A promising development covering about 5 orders 
of magnitude in time is actually taking place using the optical Kerr effect in terms of a 
pump-and-probe experiment ||129| , |13CI|| . 



Computer simulation yield as primary output. In recent Monte-Carlo | |131 
and molecular dynamics [|132| , |133|| studies the time range was extended to such long 
times, that they now deal with the same windows as achieved by the above mentioned 
modern experiments. In molecular dynamics simulations many different observables 
can be studied which are not accessible by any other experimental methods. The data 
quality is so good that the results can be compared with more sophisticated mode- 
coupling models and with numerical solutions of the mode-coupling equations, where 
the static structure factor is the only input. However, simplifled models have to be used 
approximating the real molecule ||35|, |134| , |135| , |136| , |137| , |138|| . The need for extending the 
simulations to times long enough to follow the relaxation processes in the supercooled 
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and glassy state still necessitates a compromise with respect to the complexity of the 
model and the number of interaction sites. Lewis and Wahnstrom used a three-site 
model mimicking OTP, where each phenyl ring was treated as a single interaction site 
134| , |135| , |137|| . Although simple, this model essentially reproduces the dynamics of 



OTP, in particular the anomalous behaviour of the Debye- Waller (MoBbauer-Lamb) 
factor and the existence of a fast motional process. Later, Scortino et al. extended the 
calculations of Lewis and Wahnstrom by more than one decade towards longer times 



138|| . In addition to the translational correlation functions the rotational ones were 
calculated and compared with theory. However, the planar nature of model makes it 
easier for the molecules to slide past one another. The 18-site, three-ring model of 



Kudchadkar and Wiest ||136|| and later of Mossa et al |^ takes into account the three- 
dimensional nature of the OTP molecule and may provide a better representation of the 
reorient at ional dynamics. In addition it retains some internal torsional motions of the 
lateral rings with respect to the central one. However, with the extended complexity 
only a few number of molecules could be simulated and the time range is rather limited. 
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5. Mode— Coupling Theory 
5.1. Aim 

The dynamics of normal simple liquids is governed mainly by excluded volume and 
binary collisions between molecules which, at high temperatures or low density become 
more and more statistically independent from each other. It can be described by 
generalizing the kinetic equation of Boltzmann. When the liquid is densified it is clear 
that this description breaks down at some point, as the binary collision events become 
highly correlated. This gives rise to two well known phenomena: the back-flow and 
the cage effect. Suppose a tagged-particle is trapped in a cage formed by its nearest 
neighbours which are caged themselves. On a short time scale the only possible motion 
is a rattling in this cage. The particle cannot escape and move over distances larger than 
the next neighbour distances. These correlated collisions in the cage lead to a strong 
damping of the motion and correspond to an increasing friction. On longer time scales 
and for transport over longer distances the cage must open which requires a cooperative 
rearrangement of many particles. For large densities the cage forming particles are 
subject to the same dynamics as the particle in the cage, i. e. all fluctuations have to 
be treated on the same level. 

It is the essence of the mode-coupling theory of the glass transition (MOT) that 
it indeed deals with the two phenomena cited above. Its advent has highly stimulated 
many discussions and experiments. MOT describes the evolution of structural relaxation 
in simple liquids. The cage effect is responsible for the slowing down of the dynamics 
of supercooled liquids and drives the system towards an ergodic-nonergodic transition 
at some critical temperature Tc > Tg or some critical density Uc < Ug. A smooth cross- 
over from vibrations to structural relaxation is predicted as a new dynamic feature - the 
/9-relaxation. At some point the back-flow effect counterbalances the cage effect and 
allows the system to stay ergodic until it drops finally out of equilibrium at the caloric 
glass transition. 

In the following we briefly recall the analytical solutions of the mode-coupling 
theory as derived for the idealized far as necessary for the interpretation of the 

experimental data in sections ^ and |^ It is expected that the essential features of the 



idealized theory above (or below ric) prevail in real systems ||139|| . While the analytical 
solutions for the idealized theory provide handy fit formulars, numerical solutions of the 
mode-coupling theory can only hardly be used to analyze the experimental data. The 
most adequate for our purpose is the formulation of the results in the time. 



5.2. Equations of Motion 

Starting point of the original (referred to as simplified or idealized) version of MOT 
is a generalized Langevin equation of motion for the normalized density fluctuation 
correlator $Q(t) = {pQ{0)p-Q{t))/ {pq{0)p-q{0)), where pQ{t) is the Fourier component 
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of the density fluctuation: 

d^<l>Q{t) + UQdt<l>Q{t) + Ql<l>Q{t) + Ql f dt'rriQit - t')dt,^Q{t') = 0. (14) 

Jo 

The set of integro-differential equations can be closed using Kawasakis factorization 
approximation. The kernel mQ{t) is expressed in terms of products of correlation 
functions (called the mode-coupling approximation): 

mQ{t)= Yl ViQ;h,h)^k,mk,it)=J'Q{V,^k{t)). (15) 

ki+k2=Q 

The internal control parameters V representing the strength of the coupling of different 
density modes are controlled by the external control parameters like temperature T or 
density n via the static structure factor S{Q;n,T). They can be calculated from the 
interaction potentials and vary smoothly with wave number Q and all external control 
parameters. Changing the temperature or the pressure, a dynamic phase transition from 
a liquid to a glass occurs at a critical value Vc called the glass transition singularity. 
In the vicinity of a particular value Vc a characteristic slowing down of $Q(t) is found 
0, |[. As it is not accompanied by marked structural changes we deal with a purely 
dynamic transition. 

In most real materials the knowledge of the interaction potentials is too inaccurate 
to compute the MCT solutions explicitly. However, in the vicinity of the glass transition 
singularity Vc a number of generic features of the solutions exist and are independent 
of the detailed form of the coupling constants. Most of the predictions of MCT are 
obtained from asymptotic expansions in the vicinity of the singularity. 



5.3. Solutions of Mode-Coupling Equations 
(i) The Square Root Singularity 

Denoting the distance from the singularity Vc by the separation parameter a , it follows for 
the Debye- Waller factor or the nonergodicity parameter fq = lim^^oo ^qit) in leading 
order |lil| 

. / N .c , / C'M : {a <0,ergodic) 

jQy^) ^ ]^ + Oia) : ia>0,nonergodic) ^ ' 

hq is called the critical amplitude. Let x be one external control parameter like T or n, 
then for the glass transition we have a oc (x — x^jxc- This square root cusp in fqix) is 
an important signature of Xc in the MCT. 

The Debye Waller factor can easily be measured by elastic neutron scattering. In 
practice, due to the finite instrumental resolution 8E it is impossible to measure strictly 
elastic scattering = 0. All scattering events within bE appear elastic. Therefore, 
an effective nonergodicity parameter fq is used which corresponds to the replacement 
of fqbiuS) by the area under the a peak or equivalently by the plateau height of the 
correlation function $Q(t — h/5E). 
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(ii) The [3-Relaxation 

Let us now turn to the dynamic behaviour of the density correlation function in the 

vicinity of the dynamic phase transition. The existence of a /^-relaxation with a critical 
spectrum is probably the most important result of MCT. It concerns the mesoscopic 
time range between the microscopic short time dynamics determined by VLq and the a- 
relaxation dynamics at long times. This time range corresponds to the frequency range 
around the minimum in the susceptibility spectrum ^("{uj) separating the microscopic 
peak from the low frequency a-relaxation peak. In the ideal version of MCT the shape 
of $Q(t) in this region is completely determined by a separation parameter a and a 
material dependent constant A. 

Strong predictions are made in the /3-relaxation regime, the external control 
parameters being sufficiently close to the critical point. For to = <^ t <^ t and for 
|a; — a;c| sufficiently small the correlators $Q(t) stay close to the plateau value /q. The 
equation of motion (0) can then be asymptotically developed in terms of |$Q(t) — fgl- 
Note, thereby one implicitly defines a time interval, where $Q(t) is close to fg called 
the /3-regime. One obtains in leading order ||141| , |142| 



'i>Q{t) = f^ + HQg±{t/t^), (a^O). (17) 

where the amplitude and the cross-over time 

Hq = hQCo\a\^/^ and t^ = to\a\^/^'' (18) 

depend critically on the distance from the transition point, cq is a material dependent 
dimensionless constant describing the proportionality between the theoretical and the 
experimental control parameters. 

It is a remarkable property that ^qit) factorizes into a purely Q-dependent 
amplitude Hq and a purely time-dependent function g±{t/ta). Applied to the particle 
density n(r, t) this means that the variations of the density in space are uncorrelated 
with those in time. 

The scaling function g±{t) depends only on the line shape parameter A and can be 
calculated from the expansion [|143| 

~ ^ / ± Air + ^2?" ± . . . : (to < t < t., a ^ 0) 

W ^ B,/{Bi^) + ... ■ {t^<^t<^T,a<0). 

The exponents a and b as well as the expansion coefficients Ai and Bi and the cross-over 
time are tabulated functions of A. Starting from short times t t^, g± slowly decays 
to zero. For times t longer than t^j it changes its sign for cr < and initiates the further 
decay of the correlations. 

The two non-universal critical exponents a and b are related to the line shape 
parameter A via the transcendental equation 

r^(i-«) r^(i + ^ o<a<t<i (20) 

r(l-2a) r(l + 26)' ' ' ^ ' 

from which follows that 1 /2 < A < 1 . 
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On leaving the plateau region, we obtain from (|1^) the asymptotic power law 

= fQ- BhQ\a\'/\t/t„f = fQ- BhQ{t/rf (21) 

with a new universal time scale 

r = tokr^ l=^ + h (22) 
la lb 

which leads up to the final decay of correlations in the a process. Its fractal time 
dependence is the key to stretching ||140|| . Eq. (|2T|) is called the von Schweidler law and 
its region of validity is often rather limited. However, it gives a theoretical explanation 
of the often observed time-temperature-pressure principle 

$Q(t;X) = $Q(t/r(X)), X = T,P,.... (23) 

This means in general, that the line shape may still depend on Q but not on X. A 
change of the external control parameter only causes a rescaling in time. 

MCT predicts that the relaxation times r of all correlators that couple to density 
fiuctuations should show a divergent behaviour near the critical point in the form of 
a power law with the same universal exponent 7. This property is called the a scale 
universality. The exponent 7 can be calculated once the line shape parameter A is know 
and vice versa. 

The final decay of the correlations is often parameterized by the Kohlrausch 
stretched exponential ([l|). Note, that the Kohlrausch law is not a general solution 
of the MCT equations. However, it is an asymptotic result in the limit of large wave 
numbers ||144| . 



(iii) Q-dependence 

In the framework of MCT the Q dependence of fg, hg, Aq,tq and jSq is determined 
solely by the static properties of the system and have so far been determined from 



numerical solutions of mode-coupling equations for systems like hard spheres R 145 



Lennard- Jones ||146|| spheres, and binary mixtures of hard spheres ||147|| and Lennard- 



Jones spheres |148 



In the incoherent case, i. e. the tagged-particle correlations, fq refiects the single 
particle oscillations and it closely corresponds to a Gaussian distribution in space 
fq = exp{—Q^r^) |149|| . The half width of fq is inversely proportional to the root- 
mean-square displacement of a particle. The amplitude hq starts to increase with a 
behaviour that is also expected from harmonic theory of vibrational excitations in the 
high frequency low temperature limit. 

For the coherent case, i. e. the density correlations, the amphtudes fq and hq 
refiect the oscillatory character of S{Q) and are found to oscillate in phase and out of 
phase with the static structure factor S{Q), respectively. 

Theory does not claim that the short time expansion of (|lD matches exactly (|2ll) . 
However, it makes it plausible that the a-relaxation time Tq follows the Q dependence 
of fq/hq. The a-relaxation time of the density correlations Tq is found to vary strongly 
around Qo, the position of the first structure factor maximum ||145|| . Such oscillations 
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of tq are in a more general context known as de Gennes narrowing ||113|| . Since in general 
[3q ^ h oscillations in fq and h'^ should also lead to oscillations in [3q. 

The a-relaxation time of the tagged-particle correlations Tq decreases 
monotonously close to a law while (3q varies monotonically from 1 in the Q — > 



limit to a constant value in the large Q limit ||145| . 

For both, tagged-particle and density correlations the lineshape parameter A and 
the crossover time t„ are identical and ^^independent. 

5.4- Extentions of Mode-Coupling Theory 

In the idealized theory below the a-relaxation time r and the cross-over time t^j 
are predicted to diverge at with the viscosity rj. However, in reality they stay 
finite at this temperature. The divergence is due to the overestimation of the cage 
effect in the approximation used for the memory kernel rnQ{t) (|T3|). In the extended 
theory a-relaxation and cross-over evolve smoothly with decreasing temperature below 
Tc. Including the coupling to density currents, i. e. products of type $Q(t)$Q(t) or 
$Q(t)$Q(t) in (P!3|), lifts the divergence at Tc and the system remains ergodic even 
below the critical temperature. This so called "phonon assisted hopping" processes 
counterbalance the cage effect and restore ergodicity at Tc. Thus, the ideal glass 
transition is replaced by a change of mechanism of motion from liquid-like diffusional 
motion above Tc to solid-like jump processes below Tc. 

The equations and solutions of the extended version are rather involved and we 
expect that the essential features of the idealized version concerning the /3-relaxation 
prevail. For several glassy systems, in particular OTP, we will see that nature comes 
very close to the ideal theory. 

Equations (pJ]) are the predictions of the theory about the leading asymptotic 
behaviour for the time and the temperature (or density) dependence of a generic 
correlator. Recently the next leading order corrections to the asymptotic solutions 
have been calculated ||149| , |150|| . These corrections lead to some modifications in the 



early /3-regime for t^ <^t <^ta, namely 

^Q{t) = fQ + hgito/tni + [Kq + A](toA)n (24) 

with a Q-independent constant A and a Q-dependent constant Kq; both are 
temperature independent. In the late /5-regime the correlation function is now predicted 
to behave like 

^Qit) = f'Q- hgit/rni - [Kq + A'](t/r)^}. (25) 

The corrections are the second terms in the curly brackets in (0) and (^5)). The 
important result is that the Q dependent correction constant Kq is the same in both 
equations. In addition, the constanst A and A' are related with each other. Indications 
of the existance and the correctness of these corrections have been found in molecular 



dynamics simulations [151 . 
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The success of mode-coupling theory for simple liquids has stimulated a 
considerable amount of theoretical work to extend it to molecular liquids. Recent 
contributions include the extention to the rotational dynamics of a linear molecule 
immersed into a system of spherical particles ||152|| and the generalization to molecular 
liquids of linear and rigid molecules ||153|| . Using the projection operator formalism 
an equation of motion is derived for the correlators Sim,i'm'{Q,t) of the tensorial one 
particle density pim{Q,t), which contains also the orientational degrees of freedom for 
/ > 0. Now, the translational and the orientational degrees of freedom freeze in in 
a different manner and at different critical values and demonstrate a hierarchy for the 
freezing: the orientational degrees of freedom can never freeze in before the translational 
degrees of freedom are frozen. 

MCT has also been extended to treat the dynamics of full molecular systems using 
a site-site representation ||154|| . A further step has been undertaken to describe the 
orientational and translational dynamics of liquids composed by rigid molecules of 
arbitrary shape ||155 . 



As in neutron scattering we are unable to unravel the reorient ational correlation 
functions we do not go into detail here. The interested reader is referred to the cited 
references. 

Generally, the mode-coupling approximations ( jTS]) are difficult to control as there 
is no obvious smallness parameter. Thus, MCT predictions have to be checked by 
comparison with experimental and numerical results. This will be done in the next 
sections. 
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6. Temperature Dependent Experiments 
6.1. Experimental Considerations 

As already mentioned in section ^ the crystallization tendency of OTP requires some 
care. Our first experiments have been performed with aluminum containers where the 
temperature range between 245 K and 290 K was inaccessible due to crystallization. 
Later, the crystaUization tendency could be significantly reduced using high purity 
samples in clean glass capillaries with very smooth container surfaces. For incoherent 
scattering experiments with OTP-c/q several hundred soda lime glass capillaries of 
0.32 mm inner diameter and 0.04mm wall thickness were used. For the coherent 
scattering experiments with fully deuterated 0TP-_Di4 below soda lime glass 
capillaries with an inner diameter of 1.2 mm and a wall thickness of only 10 yum were 
used. The sealed capillaries were arranged on a circle of an appropriate sample holder 
whose diameter was adapted to the neutron beam size mimicking a hollow cylinder 
geometry. This geometry keeps multiple scattering and self-shielding effects relatively 
isotropic. Furthermore, the use of capillaries allows to check visually for any traces of 
crystallization. Another positive side effect not considered before: The glass tubes can 
be used for scattering experiments beyond Q > 2.7 A^^ as no Bragg reflection occurs 
contrary to Al sample holders. 

High purity samples were obtained by repeated recrystallization out of a hot 
methanol solution and subsequent vacuum distillation. Before filling the capillaries were 
flushed with distilled water to remove dust and baked out at about 850 K for several 
hours. After filling, the samples were tempered at about 350 K for several days. When 
left undisturbed at room temperature no crystallization occured for more than 2 years. 

The measured transmission of a coUimated beam of the samples was always higher 
than 90 % which is generally regarded as a good compromise between the conflicting 
requirements of high single scattering and low multiple scattering. Therefore, no 
attempt has been made to correct for multiple scattering leading to some reservations 
in particular at low wave numbers, say Q < 0.5 A^^. 



6.2. Static Structure Factor 

Although OTP is a relatively simple molecule the static structure factor S{Q) is rather 



involved [|156|| . A straightforward connection of the structure factor peaks to single 
molecular units is no trivial task since the static structure factor S{Q) is a weighted 
sum of atomic correlations. S{Q) of supercooled deuterated OTP as measured on the 
diffractometer D20 of the Institute Laue Langevin (ILL), Grenoble, is shown in figure 0. 
In contrast to atomic systems the main peak of the static structure factor is splitted into 
two maxima at about 1.4 and 1.9 A~^. With decreasing temperature a slight shift of 
the peak positions to higher Q values due to the increase in density n and in addition an 
increase of the peak height at 1.9 is observed. The peak height at 1.4 is slightly 
reduced which is probably connected to the decrease of the isothermal compressibility 
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Figure 13. (a) Static structure factor of fully deuterated OTP at 314 K (full line) and 
255 K (dashed dotted line) as measured on the neutron diffractometer D20 of the ILL, 
reproduced from |156|. (b) Static structure factor in the low Q region of deuterated 
OTP from IN5. Note the appearance of a clear shoulder around Q — 0.85 The inset 
shows the temperature dependence of S{Q) at selected Q values which are indicated by 
the arrows in the figure. S{Q) shows strong temperature variations at Q = 0.85 
and Q = 1.95 A-\ Taken from El. 



nxTkBT = S{Q = 0). 

Due to experimental limitations the Q region below about 1 A^^ was not accessible 
in the D20 experiment. Therefore, we show in the right part of figure [T^ the static 
structure factor S{Q) in the low Q region as measured on the time-of-fiight spectrometer 
INS of the ILL The double peak structure is resolved and the temperature 

dependence of S{Q) agrees well with the one measured on D20. One recognizes a clear 
shoulder around 0.85 with an increasing peak height with increasing temperature. 
In inelastic X-ray scattering experiments a real prepeak at Q ~ 0.85 A^^ was observed 
||157|| together with a main peak at Q ^ 1.45 A~^ while nearly nothing is found around 
Q ~ 1.9 A-^ Note that in neutron scattering on fully deuterated OTP both Deuteron 
and Carbon atoms are seen with about the same scattering cross sections (cf. table |l]), 
while in X-ray scattering only the distribution of Carbon atoms has to be taken into 
account. 

In simulations simplified models are used for the complex molecule to keep the 
computing time reasonable. Lewis and Wahnstrom employed a three site complex, 
each site playing the role of a whole benzene ring ||135|| . The static structure factor 
of the center of mass shows a maximum around 0.8 ... 1 A^^ while the one of the sites 
is solely peaked at 1.5 A~^. A refinement of this model with 18 sites on 3 rings was 
later studied by Kudchadkar [p.36|| . The static structure factor shown in figure |T^ now 
well reproduces the experimental one, in particular at higher Q including the hump at 
3A~^ and the prominent peak at 6A~^. For smaller Q it yields, besides a double peak 
structure at 1.4 A~^ and 1.5 A~^, also a clear shallow peak around 0.8 A~^ as illustrated 
in the inset of figure |l^. The fact that the second peak of the double peak is at lower 
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Q than observed experimentally is probably due to the fact of "fusing" the carbon and 
the attached hydrogen atoms into one "united atom". In a recent work on a flexible 
molecule model a good agreement among the molecular dynamics structure factor 
and that measured by neutron scattering was obtained in the high-Q region. A clear 
shoulder was observed at 0.9 A~^, a diffraction peak at 1.5 but no double peak 
structure, only a small bump around 1.8 A^^. Their analysis in terms of a calculated 
structure factor on the molecular and ring center of mass clearly support our simplified 
picture that follows below. 
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Figure 14. Calculated static structure factor for the 18-site model of OTP. The inset 
displays the small Q region in more detail. Note the shoulder around Q — 0.8 A^^. 



Reproduced from |13£] 



In this context we point out that the structure factors of fully deuterated benzene 
and hexafiuorobenzene, that can be considered as unconnected subunits of OTP are 
nearly indistinguishable from the one of OTP both showing the double peak at the 
right positions as well as the hump at SA^^ and 6A^^ ||158||. Similar double-peak 



structures are found even in similarly simple phenyl systems like m-fiuoroanilin |159 



In this context we point out that the structure factors of fully deuterated benzene 
and hexafiuorobenzene, that can be considered as unconnected subunits of OTP are 
nearly indistinguishable from the one of OTP both showing the double peak at the 
right positions as well as the hump at 3A~^ and 6A~^ [|158|| . Similar double-peak 



structures are found even in similarly simple phenyl systems like m-fluoroanilin [|159|| . 
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The following simplified picture may therefore be developed: The peak at Q2 — 
1.9 is probably built up mainly by correlations within the phenyl rings, while the 
maximum at Qi ^ 1.45 might be associated with correlations between phenyl rings. 
The shoulder around Qo — 0.85 is presumably the most direct manifestation of 
intermolecular correlations I^S], |135| , |136|| . Its position is indeed compatible with the 
inverse van der Waals radius of rw = 3.7A[^. We don't believe that the shoulder in 
OTP is due to pure orientational correlations that can lead to a prepeak at low Q [|16CI|] . 



6.3. Square Root Singularity 

The first experimental indications of an anomalous behaviour of the Debye- Waller factor 
fq = exp{—2W) have been observed long before MCT of the glass transition was 
developed in MoBbauer-spectra of ^^Fe in ferrocene dissolved in butylphthalate ||161 



and OTP | 162 |. The results for OTP at Q = 7.3 A^^ are shown in figure One 
notices a linear temperature dependence at low temperatures. Above Tg a rather strong 
increase of 2W is observed. 
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Figure 15. Exponent of the Debye-Waller factor of ^^Fe in ferrocene dissolved in 
OTP as a function of temperature. Based on values given in |162| . The straight line 
is a linear fit to the lowest two values and zero to give a consistent picture with the 
results discussed in section |8[ Note a nearly linear temperature dependence at low 
temperatures followed by a strong anomalous decrease above Tg. 



The predicted square root singularity (|lq) was extensively tested by neutron 

IMI, |l6l, |l6l, |166 



scattering experiments |42 



At low temperatures /g can be 
easily measured by the elastic intensity in a BS experiment or by determining the plateau 
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of the vs. \ogt curve. In figure |T6| we give results for temperature dependent 

incoherent and coherent scattering. For low temperatures we simply follow the harmonic 
evolution of the mean square displacement (In /q(T) oc T). The anomalous decrease 
connected with the glass transition starts somewhere around Tg in the glassy phase due 
to anharmonic motion on a microscopic length scale. It is clear that the structural 
relaxation which at Tg is of the order of seconds cannot be responsible for this effect 
and therefore Tg as conventionally defined has nothing to do with it. 
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Figure 16. (a (*)) Incoherent Mofibauer-Lamb factor of protonatcd OTP for several 
Q- values | |163{ . A harmonic behaviour at low temperature is followed by the onset of 
an anomalous decrease around Tg. Note the cusps in the temperature dependence at 
Tc- A square-root law fit gives a common Tc ~ 290 K for all Q. (b (*)) Temperature 
dependence of the Debye- Waller factor of the density correlations for various values of 
Q. Lines are fits with a square root law leading to a critical temperature of Tc ~ 290 K. 
Towards higher temperatures a weak linear temperature dependence is assumed. The 
inset demonstrates the square root law in linearised form. 



On approaching Tc the sharp transition predicted by the idealised version of the 
theory is smeared out through thermally activated hopping processes, which are present 
in any molecular liquid. Clearly, between the caloric glass transition temperature Tg 
and Tc structural relaxation is present. When its time scale falls into the experimental 
window quasielastic broadening appears and its total intensity can no longer be 
determined from strict elastic scattering alone. Instead, the integral over the whole 
a peak has to be taken. The results are presented in figure |1^ showing all features 
predicted by ([16|). 

An important question to answer is whether the anomaly of the Debye- Waller factor 
may be caused by any intramolecular motion such as the librational dynamics of the 
lateral phenyl rings relative to the central one. This can be checked by measurements 
of the Debye-Waller factor in isotopically substituted OTP systems, namely fully 
protonated OTP-do, fully deuterated 0TP-(ii4 and selectively deuterated at the lateral 
phenyl rings OTP-^iq. Assume two types of motion, a center-of-mass (COM) motion 
and an intramolecular (IM) motion of the lateral phenyl rings. Then the central ring 
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only contributes to the COM motion, while the lateral rings contribute to both the 
COM and the IM motion. From the total scattering cross sections (cf. table |I]) for the 
nuclei participating in each motion it follows, that in particular in OTP-rfio neutron 
scattering is nearly insensitive to the lateral rings (they are "hidden") as the dominat 
scattering steems from the central ring protons. Consequently, if the /3-process, which 
is responsible for the anomalous decrease of /q(T) is due to intramolecular motion, 
it should be much less pronounced in OTP-diQ. As the Q dependence of the DWF 
is essentially the same for the three isotopically substituted OTPs [[101|1 and can be 
approximated by a Gaussian fq = exp[—^Q'^{r'^)] we can deduce the mean square 
displacement shown in figure |l^ From the comparison of the observed anomaly in 
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Figure 17. Mean square displacement of three isotopically substituted OTPs obtained 
from elastic scans on IN13. Within experimental error the three data sets are identical. 



Taken from [101 



OTP-do and OTP-ciio, we can exclude any intramolecular phenyl ring dynamics as the 
dominant mechanism for the f3 process. The data thus support MCT proposing a center- 
of-mass motion as precursor of the glass transition. Note, that although the values for 
(r^) of 0TP-(ii4 conicide with those of OTP-do and OTP-dio, the comparison is biased 
as in 0TP-(ii4, which is a dominant coherent scatterer, we do not observe the single 
particle mean square displacement. In addition, with a higher Q-resolution one can 
clearly observe oscillations in fq (see figure ^2|) leading to deviations from the simple 
Gaussian behaviour. 

The anomalous behaviour and cusp of the nonergodicity parameter fq has also been 



observed in computer simulations [|134|, |135|] . The cited critical temperature = 280 K 



is in reasonable agreement with the one determined from neutron scattering. X-ray 
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results allowed for the identification of a temperature — 290 K and a cusp from the 
temperature dependence of a high energy propagating soundlike mode ||120|. 



6.4- Critical Correlations 

On BS spectrometers only indirect information can be gained about the existance of a 
fast process. In order to study the (3 process directly shorter times have to be accessed 
which is possible on TOF spectrometers. Whether data are analyzed in frequency or 
time is a question of convenience. 

One remarkable prediction of MCT is the factorization property of the /3-relaxation. 
In a certain frequency and temperature or pressure range, all observables are expected 
to have the same spectral distribution. This implies that the dynamic structure factor 
S{Q,Lj) factorises into a Q- and a u dependent part as a direct consequence of the cage 
effect. Figure |T^ demonstrates the validity of the factorization property for incoherent 
and coherent scattering: the rescaled intensities S{Q,u!)/ S{Q)Aq = G^u) fall onto a 
common master curve for different wave numbers Q p2| |107| , |167| | . 
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Figure 18. Test of the factorization property for coherent at 293 K (a) and 
incoherent scattering data at 324K and 60MPa (b (*)). The rescaled intensities 
S{Q,uj)/ S{Q)Aq = G{uj) for different wave numbers ranging from 0.4 to 2A-^ 
are shown to fall onto a common curve in the energy range between 0.1 meV and 
1 meV. For comparison a power law to^^'^, a — 0.3 is shown which is the proposed line 
shape from theory. The steeper lines are the measured resolution functions and set the 
lower bounds. The insets show the scaling factors Aq which are proportional to the (3- 
relaxation amplitude hg. Note for coherent scattering the two minima corresponding 
to the prepeak and the first maximum of the static structure factor. 



For a detailed study of the /5-relaxation we Fourier deconvolute the dynamic 
structure factors as the quasielastic broadening is of the order of the instrument's 
resolution. For an analysis in time domain we combine data from different spectrometers 
with overlapping time ranges to enlarge the dynamic range. Some results for the tagged- 
particle ||107|| and the density correlations |164|| at different temperatures are shown 
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in figure In Ref. ||167|| a second Fourier transform was liereafter performed to analyze 
tfie scaling behaviour of the dynamic suceptibility x"{^) around its minimum. 
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Figure 19. Normalized intermediate scattering fmictions ^{Q,t) for density (a (*)) 
and tagged particle fluctuations (b) (taken from [^) for two different wave vectors at 
various temperatures. The ranges of the instruments used are indicated. Solid lines 
are fits with the MCT scaling law (p7|). The line shape parameter A = 0.77 was kept 
fixed for the incoherent data while for the coherent data is was fitted freely. The decay 
of the correlations is described over several decades. For the coherent data we included 
KWW fits (0) to the long time tail of the correlation functions (dotted lines). 



For fits with (P!7| ) we use the tabulated expansion coefficients ||143|| of the scaling 
function g\{t/ta) which is completely determined by the exponent parameter A. 
Meaningful four-parameter fits are often not warranted by the dynamic range covered, 
here in particular for the incoherent data shown in figure |19| on the right. However, 
as the line shape is not very sensitive to small variation of A it was kept fixed at the 
value 0.77 previously derived from a scaling analysis of the long time asymptote of the 
correlation function <^{Q,t) ||163| , |168|| . This value is compatible with the one from the 
power-law fit to the viscosity ri{T)/T in figure ^ ||168|| . For the coherent data by the 
combination of time-of-fiight (INS), spin-echo (INll) and back-scattering (IN13, IN16) 
data a time range of now more than four decades is covered as can be seen in the left 



part of figure |T9|. INS, INll and IN13 have a broad overlap and their $((5, t) show very 
good agreement after the absolute scales were adjusted by typically S%. IN16 data 
couldn't be obtained in absolute units and were matched to INll or IN13 data. Note, a 
combination of spin-echo data with TOF and BS data has rarely been attempted 
possibly due to the different weighting of coherent and incoherent scattering in (p!2|). 
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However, when the coherent scattering largly dominates the signal, for example around 
the static structure factor maximum such a combination is reasonably possible. Now a 
four-parameter fit seems more promising and the results for the line shape parameter 
A is shown in figure |20|. The values for different temperatures T and wave numbers 
Q indeed scatter statistically around the mean value A = 0.77 validating our previous 
analysis P|, [TU71. 
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Figure 20. (a (*)) Line sliape parameter A from fits to the density correlation 
functions (figure ^9|) which have been obtained by combining up to four different 
instruments (IN5,fNf f ,IN13,rNf 6). The values lie between 0.7 and 0.8 and scatter 
statistically around a mean value of 0.77 which is consistent with the one determined 
by other means. Note, no significant temperature (different symbols) or Q^dependence 
is observed in accordance with the prediction of MCT. (b (*)) Von Schweidler exponent 
b for density correlations from fits to the master curves. The data scatter statisctically 
around a mean of 0.53 which is consistent with A = 0.77. 



The /9-relaxation scaling is expected to hold only in a wave number and temperature 
dependent time window [tmm{Q,T),tjnax{Q,T)] around the cross-over time to-. As the 
scaling law is a long-time expansion where details of the vibrational motion are neglected 
it should fail to hold for times shorter as or of the order of a typical vibrational period, 
say a psec. Towards longer times the a-relaxation which itself strongly depends on Q 
and T starts to merge with the late f3 regime and renders the identification of a fit 
range rather ambiguous. As the predictions of MCT are based on expansions around 
we expect the scaling law to hold best close to the critical temperature. There, on the 
other hand, additional transport processes not included in MCT's elementary version 
may obscure it. 

From figure we can see that the fits with ( [T7| ) are able to decribe the decay of 
tagged-particle as well as density correlations over several decades in time. It might 
come as a suprise that the correlators at the highest temperatures are described down 
to about /q/2 as the leading order solutions ^{Q,t) ( |17|) have been obtained by an 
expansion in the distance from the plateau \^{Q,t) — 

The physical meaning of such fits has to be judged from the comparison of the 
obtained parameters with the theoretical predictions (0). Such comparisons should 
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comprise not only temperature T and / or pressure P but also the wave number 
Q dependence. The temperature dependence of hq and t^ for tagged-particle and 
density correlations is show in figure ^ confirming the asymptotic prediction ([T8|). 
Extrapolation of the T-dependence gives consistently Tc ~ 290 K for all data sets 
coinciding with the critical temperature determined from the square root anomaly of 
the Debye- Waller factors and the divergence of the viscosity. 
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Figure 21. Temperature dependence of the amplitude Hq and the cross-over time 
as obtained from fits to the correlation functions ^{Q, t). The parameters are plotted 
such that a linearized T dependence is expected according to (p^. Coherent (a,b) as 
well as incoherent (c,d) data extrapolate consistently to Tc ~ 290 K. While the absolute 
values of both, the incoherent and the coherent cross-over times t,, are quite the same 
the amplitudes Hq differ. Taken from 107 1. 



Also the wave number dependence shown in figure at least qualitatively 
agrees with results from numerical solutions of mode-coupling equations f^, |145|| . For 
incoherent scattering /g reflects the single particle oscillations and it corresponds closely 



to a Gaussian distribution ||163| , |168|| . For coherent scattering the plateau height /g and 
the amplitude hq are obtained from a fit by (|16D to the DWF shown in figure |16| and from 
a fit of ( ]T7| ) to t) shown in figure |T9|. Both results are nearly identical and oscillate 
in and out of phase with the static structure factor S{Q), respectively. Remember, the 



scaling factors Aq in the insets of figure ^ are proportional to hq, which independently 
gives strong support for our findings. 

When the abscissa of the theoretical curves [^, |169|| is scaled to a hard sphere 
diameter of a ~ 7.5 A even a more quantitative agreement is found for both, the coherent 
and incoherent results. 
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Figure 22. (a (*)) Wave number dependence of the non-ergodicity parameter fg 
(open symbols) and the /3-relaxation amplitude hg (closed symbols) for the density 
fluctuations compared to the static structure factor (line). The results stem from the 
fits to the correlation function (squares) and to the DWF (circles), /g varies in phase, 
hq out of phase with the static structure factor, (b (*)) Same as in (a) for the tagged 
particle fluctuations. The lines are predictions of theory redrawn from |^ with the 
abscissa scaled to a hard sphere diameter a ~ 7.5 A which corresponds to the prepeak 
of the static structure factor S{Q). 



The factorization property has also been verified for the coherent intermediate 



scattering function in an appropriate time regime [42|. According to (O) the quantity 



($(Q, t) — fg) I hq is independent of Q. This has been tested using the fit parameters fq 
and hq for wave numbers between 0.8 and 2A~^. The cross-over time t^ and the 
time range over which the factorization is fulfilled decrease with increasing temperature. 
This factorization is only possible if the line-shape parameter A and the cross-over time 
to- of the /5-correlator are both independent of Q. A fit of g\{t/t^) to the combined data 
gives A = 0.78 in excellent accord with the values determined by other means [Q. 

Here we want to follow another approach. As already mentioned in section ^ in the 
context of equation (p!7D, the whole time dependence of any arbitrary time-correlation 
function ^q{t) which couples to density fluctuations is given by the Q-independent 
function g_{t/t„){a < 0). Introducing the following function ||151j| 

- - (^"^ 

one can test this prediction even without fitting with (|13). Here t' ^ t" are abitrary times 
in the /3-relaxation regime. From ([T7| ) we immediately recongnize that the function 
Rq{t) is independent of the correlator, i. e. of Q. To test this property we take the 
density correlation function at 293 K of figure For the five different Q (0.85 A^^ < 
1.8 A^^) we determined Rqif) choosing arbitrarily t' = 5ps and t" = 50 ps. Figure 23 



shows that in the /5-relaxation regime the different Rq{t) indeed collapse onto a master 
curve demonstrating the validity of the factorization property. This is not a trivial result 
as outside the /3-relaxation regime the curves fan out. One can even try to go further. 
Regarding the next order corrections (|2^) and (|25|) it can easily be shown that for the 
function Rq{t) the order of the corrections in the early /9-relaxation regime and in the 




Figure 23. (*) Time dependence of the ratio RQ{t) given in ([2^ ) demonstrating the 
factorisation property for T=293K. In the upper right inset Rgit) for fitted curves is 
shown. For the lower left inset see text. 



late /3-relaxation regime are related in the following way: Counting from top to bottom 
the first (second, . . .) curve at short times (label i) is also the first (second, . . .) at large 
times (label j) which means that i = j. Both times have to be taken slightly outside the 
range where (p!?]) holds. Counting at t ~ Ips and at t ^ 1000 ps gives us the function 



displayed in the inset of figure |26|. Although only five values are considered here 
there seems to be no clear correlation i = j but rather an anticorrelation. This is in 
contrast to results from molecular dynamic simulations p51| j. As the experimental data 
scatter quite large we fitted the correlation functions with a power law similar to the 
(P^) however without any constraints to have a good smooth representation of the data 
in the time range from about 0.2 ps to about 2000 ps. We then calculate Rqit) for the 
fitted curves which are displayed in the upper right inset of figure |2^. Identical results 
are obtained concerning the factorization and the function as from the experimental 
data. Here again the limitations of neutrons scattering show up: To tackle this problem 
properly, high energy transfers (very short times) in combination with high energy 
resolution (long times) (preferably on one spectrometer) would be needed in addition 
to a good Q resolution and a large wave number range to have a considerable amout of 
different correlation functions. 

The fast dynamics has also been experimentally observed by depolarized light 
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scattering |125| , |127|| . Comparisons with idealized MCT yields a critical temperature 



Tc ~ 288 - 292 K with line shape parameters of a = 0.33, b = 0.65, A = 0.7 [^, |T25 



and equally ~ 288 - 292 K with a = 0.32, b = 0.61, A = 0.72 [|12^. The numbers for 
Tc and A are consistent within the experimental uncertainies with the values used to 
describe the neutron scattering results. Depolarized light scattering has a rather large 
dynamic range covering about 4 decades and high statistics such that a comparison 
with the extended MCT became possible [|1271] . The analysis yields T^. ~ 276 K with 



a = 0.3,6 = 0.54, A = 0.76. 

The dynamics below the critical temperature was investigated less extensively. 
Below Tc the idealized MCT predicts a nontrivial behaviour of the generalized 
susceptibility x"{^)- ^ white-noise spectrum at very low frequencies o*^ ^ should 

cross-over to a power law x"{^) oc uj"',a < 0.4 resulting in a "knee" in x"{^) 
intermediate frequencies. Rather counterintiutive is the nontrivial prediction for the 
temperature dependence of the position of the "knee" -frequency: tu^^ should increase 
with increasing temperature and diverge at Tc. In our neutron scattering experiments 
in the temperature range 200 < T < Tc no indication for such a cross-over is observed 
in the frequency range u > 0.1 meV (25 GHz). The same observation was made in 
careful light scattering experiments ||170|| . However, this does not necessarily mean that 



the "knee" does not exist as the influence of the hopping processes i. e. its strength 
is unkown at present. If the hopping is sufficiently small the ideal spectrum is rather 
undisturbed and the "knee" should appear, however, if the hopping term is large the 
inital a decay can make it dissapear. 

So far, the only support of the existence of such a diverging time scale below 
Tc comes from PCS experiments on a hard sphere colloidal liquid [0) i|) 01 and on a 
polymer micronetwork colloidal liquid [ |171[ | where no microscopic dynamics and no 
hopping processes disturb the /3-process. 

6.5. a -Relaxation 

The a-relaxation can be directly studied e. g. either by energy resolved BS or NSE 
experiments. For instance in figure |2^ we show the coherent scattering function ^{Q,t) 
of 0TP-(ii4 as measured on the SE instrument INllC using the multidetector set-up. 
For comparison, S{Q, uj) as measured on INIOB with a heated monochromator is shown 



on the right hand side of figure |2J for several wave numbers Q. In both cases the 
line-shape cannot be described by a single Debye function but can be well fitted by a 
Kohlrausch function. 

The validity of viscosity scaling and time-temperature superposition is still an 
ongoing debate. We have seen that it is a very good first approximation at least in a 
narrow temperature range. Statistics and the limited dynamic range of one instrument 
hardly allows for a serious test. However, consistent results with acceptable error 
bars can be obtained by combination of measurements from different temperatures (or 
pressures). Assuming that viscosity scaling holds over the covered T-range we employ 
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Figure 24. (a (*)) High resolution spectra of OTP-du at Q — 1.46 A^^ for several 
temperatures as measured on the SE instrument INllC using the multidetector. The 
fit curves demonstrate that the line shape of the structural relaxation is well described 
by a KWW function according to (|l|) with /3 ~ 0.57 and not by a simple Debye process 
drawn for comparison (dashed line), (b (*)) Dynamic structure factor S{Q,uj) at 320 K 
for Q = 0.85 (□), Q ^ 1.45 A^^ (■) and Q = 1.90 A"! (0) as measured on the 
BS instrument INIOB with a heated monochromator. Lines are fits with a Fourier 
transformed Kohlrausch function with Pq = 0.57 fix for all Q. 



the following interative procedure to construct master curves ^{Q,i): We use published 
viscosities [|^, |2^, |28| to rescale experimental times to 

i = t[v{n)/n]MT)/T] (27) 

with an arbitrary normalization Tq = 290 K. From Kohlrausch fits to ^{Q,i) we obtain 
first approximations for the line shape parameter /3q and the relaxation time tq{To). By 
means of this tq(Tq) the unsealed tq(T) of the individual temperatures are calculated 
and are then kept fixed together with the line shape parameter (3q in fits to the individual 
^{Q,t). The resulting set of amplitudes Aq is used to obtain reduced correlators 
$((5, t)/AQ and the procedure is repeated. It is numerically stable and converges after a 
few iterations ||163|| . To get the temperature dependence of each parameter the individual 
^{Q,t) at each temperature can then be fitted with one free parameter while the other 
two are kept fixed at the values determined from the iterative procedure. 

The results in figure ^ confirm the validity of viscosity scaling for temperature 
dependent tagged-particle and density correlations. Note, that the scaling is a necessary 
but not sufficient condition for the applicability of MCT. 

The temperature dependence of the relaxation time tq of the density correlations at 
the structure factor maximum is shown in figure 26. Combination of data from several 
spectrometers using the iterative procedure described above makes it possible to cover 
nearly five decades in relaxation time. The structural relaxation time follows nicely the 
viscosity confirming the validity of viscosity scaling. 
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Figure 25. (a) Master curves <t{Q,t)/AQ{T) for the a-relaxation of the density 
fluctuations at different temperatures, obtained by rescaling times with, viscosity 



r]{T)/T. From [f72|. (b) Typical master curves obtained using the scaling assumption 
for the tagged particle fluctuations. From |162]. 




Figure 26. (*) Temperature dependence at ambient pressure of the a-relaxation time 
TQ at the structure factor maximum Q = 1.4 obtained from the scaling analysis of 
BS (□), NSE (0: single detector, ♦: multidetector) and TOF (■) data. The solid hue 
is the scaled viscosity interpolated from literature data |2^, ^ ^ . The inset shows 
Tg^^^ vs. T from BS and TOF data at three different Q-values using 7 = 2.59. A 
linear behavior extrapolating to zero at Tc ~ 292 K is observed for all Q in accordance 
with the prediction of MCT. 
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In the inset of figure ^ we test the divergence of the a scale predicted by MCT 
(p^). Assuming |cr| oc (T — Tc) a plot of Tq^^"' vs. T should produce a straight line 
extrapolating to zero at Tc. Indeed, a linear behaviour is observed for all wave numbers 
studied extrapolating to zero at Tc ~ 292 K. This also supports the a-scale universality. 
The deviations from this law in a narrow temperature range close to Tc may be due to 
thermally activated hopping processes not contained in MCT's ideal formulation and 
smearing out the sharp transition. As a consequence, testing the power law only in a 
narrow range around the critical point could give a wrong idea of, e. g., 7. 

The cross-over from the (3- to the a-process is described by the von Schweidler 
law (|2l|) . Thus, it should be detectable in the short time behaviour of the a-process. 
By plotting the master curves in the form log[l — ^{Q,i) / fq] vs. logt this should yield 
straight lines with slope b for short rescaled times. The result for coherent and incoherent 
data is shown in figure ^ Both data sets show a von Schweidler behaviour at short 
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Figure 27. Fits of the von Schweidler law to some linearized master curves for 
coherent scattering (a (*)) and incoherent scattering (b) (from |163|) in order to show 
the fractal time dependence t''. The slopes b scatter around a mean value of 0.53 for 
both coherent and incoherent data. This value is in excellent accord with the value 
imposed by the line shape parameter A=0.77. 
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times while at long times deviations occur. In accordance with the theory the values of 



the von Schweidler exponent b do not show any systematic Q dependence (see figure pO| ). 
They scatter around a mean value of 6 = 0.53 ± 0.03. 

In order to investigate the Q dependence of the parameters describing the density 
correlations of the a-relaxation a high Q resolution is necessary. For the question studied 
here spin-echo seems less appropriate than high resolution multidetector instruments as 
the Q resolution is implicitly limited by the transmision band of typically AA/A 10 %. 
Note, increasing the Q resolution by use of a narrower bandwidth or a graphite 
monochromator leads to an unacceptable reduction in intensity in particular if the 
experimental beam time comprises several Q and T scans. 

The big advantage of spin-echo, the large dynamic range of about three decades, 
can be compensated either by combining data from time-of-fiight and backscattering 
spectrometers or by the use of master curves. On the BS instruments IN16 (figure ^) 
and INIOB (figure |2^) the Q resolution was improved by covering part of the analyzer 
crystals by large cadmium shields ||1 72|| . This reduces considerably the number of 
detectors which can be used in the analysis. Multidetector TOF instruments intrinsically 
have a good Q resolution. However, now a compromise between wave number coverage 
and energy resolution has to be chosen. On INS the elastic energy resolution (fwhm) 
was set to 25 /ieV restricting the accessible wave number range to g < 1.7 A-\ As 
25 /ieV is a rather broad resolution for a study of the a-relaxation higher temperatures 
are required to make the dynamics appear in the experimental energy window. 
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Figure 28. Density correlation function, obtained by Fourier deconvolution of 
colicrent neutron scattering spectra measured on the time-of-flight spectrometer INS. 
(a) for different temperatures at the structure factor maximum 
Temperatures are from top to bottom: 313 K, 320 K, 330 K, 340 K, 360 K, 380 K, 
400 K. (b) on an enlarged scale for 360 K at different wave numbers Q. Solid lines 
are Kohlrausch stretched exponentials as determined from a master curve analysis. 
Without any fitting one recognizes that correlations decay steeper at Qi than at other 
wave numbers. Taken from |172|. 
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Some of the data are shown in figure p^a. Up to high temperatures, far above the 
melting point, correlations decay in two steps. The second step is strongly stretched 
and temperature dependent as expected for the a-relaxation. In figure p8| b the long 
time tail belonging to the a-relaxation is plotted on an enlarged scale for different 
wave numbers Q. Not only the relaxation time, but also the line shape exhibits a 
pronounced Q dependence. At Qi = 1.45 A~^, the position of the first maximum of 
the static structure factor, the correlations decay steeper, i. e. with less stretching, 
than at any other wave number. Such a behaviour can also be seen in the master 
curves displayed in figure |2^. As one can easily imagine from the data in figure p8| a the 
individual correlation functions $(Q,t;T) do not allow for a three parameter fit of (|lD: 
Either, for low temperatures the a-relaxation is outside the accessible time window of 
the spectrometer, or at high temperatures, where the decay is fully in the window it 
is heavily disturbed by the /3-relaxation. Then a separation of both contributions is 
not obvious and the fit range is somewhat arbitrary. Note, that even at the highest 
temperature the decay of the correlation function cannot be fully described by a single 
Kohlrausch law. Moreover, a free fit with (|lD yields unphysical results like a strongly 
decreasing stretching to values as low as 0.3 and an increase of the amplitude with 
increasing temperature. 

For a more quantitative analysis we employ the time-temperature superposition to 
obtain master curves. Indeed, without further adjustment, the data converge towards 
a temperature independent long-time asymptote ||172|| . The results from this scaling 
analysis are shown in figure The amplitude factor Aq and the line shape parameter 
(5q are quite independent of temperature in the range from 313 K to 400 K in accord 
with results from light scattering [ |127]| . In contrast to an often held belief, the stretching 
exponent does not increase with temperature towards a diffusion limit /5 = 1 in the 
normal liquid phase. Even 70 K above the melting point the correlation functions $((5, t) 
remain stretched. The structural relaxation time follows the viscosity and has a generic 
Q-dependence which confirms the a scale universality. 

All parameters describing the density correlations exhibit pronounced characteristic 
oscillations in phase with the static structure factor S{Q) (figure ^9]). The amplitude 
Aq is in almost quantitative agreement with the Debye-Waller factor /q previously 
obtained from the /3-relaxation plateau (figure |2^). Note, fq and Aq not necessarily 
have to be identical. The relaxation time tq shows two plateaus around 0.8 and 
1.4 A~^ on top of an overall decrease. More important, the stretching exponent fSq shows 
a systematic variation in phase with S{Q) being largest at the structur factor maximum. 
Even around 0.8 A^^ a slight variation is suggested from our data. At 1.4 A~^ the results 
are obvious from the raw data by plotting the Q dependent decay curves in one diagram 
(see figures. ^ and ^Sp. They are fully confirmed by fits of (|T|) to numerical solutions 
of MC equations for the hard spheres system ||145|| and a soft spheres mixture ||173| . 

Some reservation is necessary for the smallest wave numbers Q < 0.5 where 
Aq starts to tend to 1. Even for 99% deuteration the incoherent scattering amounts to 
about 15 % to the total scattering cross section. In addition, in the small Q region the 




Figure 29. Kohlrausch parameter at seven different temperatures ranging from 313 K 
to 400 K. In each of these fits, one parameter was free while the two others were kept 
fixed at the values obtained from the master curve analysis, (a) The amplitude Aq and 
(b) the stretching exponent Pq are temperature independent confirming our scaling 
approach, (c) The a-relaxation time tq follows for all Q the viscosity and exhibits a 
generic Q dependence. The data ar 330 K have been omitted for clarity as at 320 K the 
results from the BS instrument INIOB are included (a). The dotted line represents 
the quantity S{Q,S20K)/Q'^ know as de Gennes narrowing. Based on [ |172| . 



coherent scattering itself should be weak as it is roughly proportional to S{Q) which 
tends towards a constant level. Therefore, the relative incoherent contribution might be 
even much stronger than 15 %. This effect can be studied in more detail by a separation 
of coherent and incoherent scattering by polarization analysis, cf. figure |12[ From this 
experiment we learn that the coherent /q or the a-relaxation amplitude Aq must be 
smaller than the incoherent one at Q < 0.7 and remains larger for higher Q values. 
This behaviour is indeed observed in figure ^ and we can conclude that at the lowest 
Q values contributions from multiple scattering cannot be excluded. They are now even 
more complicated to treat as coherent and incoherent scattering are a weighted mix. A 
separation of coherent and incoherent scattering from very thin samples i. e. with high 
transmision seems the most straightforward access to solve this problems. However, as 
outlined above quantitative results are still out of reach within an acceptable measuring 
time with the presently available neutron flux and spectrometers. 
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In a recent work on Toluene, the structural relaxation time tq followed nearly 
perfectly Scom{Q)/Q^, where Scom{Q) is the center-of-mass static structure factor 
calculated by Monte-Carlo simulation ||174||. Such a dependence, in a more general 



context known as de Gennes narrowing ||113|| , indicates the dominance of the center-of- 



mass motion, which is of purely translational nature. We therefore included in figure ^9 
S{Q)/Q'^. Around Q = 1.4 the variation of tq with Q is indeed correclty described 
by S{Q)/Q'^. The simulated center-of-mass static structure factors for OTP is peaked 
around 0.8 A^^ [[136|1 . Therefore, the plateau in tq around 0.8 A^^ strongly supports our 
simplified picture that was drawn from S{Q) in section |6.2| . 

The wave number dependence of the tagged-particle correlations measured by 
incoherent scattering has already been discussed extensively [|163| , |168|| : all parameters 
depend monotonously on Q: PindQ) ~ const., Aq oc /q decreases like exp (— cQ^) and 
TinciQ) decreases close to a Q'"^ law which is prefectly confirmed by a recent molecular 
dynamic simulation ||35[| . 
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7. Experiments under Pressure 
7.1. Experimental Considerations 

For practical reasons most investigations on the glass transition concentrate on 
temperature effects at ambient pressure Pq = 0.1 MPa although there is a clear interest in 
studying supercooled liquid dynamics in the full T, P parameter space. We emphasize 
the difference between the effect of pressure and temperature. The primary effect of 
pressure is to change interatomic distances; changes, for example in atomic vibrations 
are a secondary effect via anharmonicity of the potential. In contrast, the primary effect 
of temperature is to alter the atomic vibrations; changes in interatomic distances are 
now a secondary effect. Pressure and temperature are thermodynamically independent 
variables even though their effects on one set of observables might be the same. Pressure 
experiments have demonstrated that the density is not the only driving force for the 
glass transition ||175| , |176| , |177| , |178|| . In the microscopic approach of MCT, both pressure 



and temperature control the dynamics via variations of the static structure factor S{Q) 
which has been measured on fully deuterated 0TP-(ii4. All other experiments have been 
performed with fully protonated OTP, thereby measuring the tagged-particle dynamics. 

While low temperature is regularly used in neutron scattering, the use of high 
pressure is less standard as experiments under pressure are much more difficult to 
perform. There are special requirements for high-pressure cells used in neutron 
scattering. For example, the flux rates of current neutron sources require large sample 
volumes in order to get acceptable signal-to-noise ratios in reasonable time. On the 
other hand the cell body in the neutron beam gives rise to unwanted scattering and 
absorption which demands a minimum of material in the beam. 

High pressure experiments are possible only in a limited temperature-pressure range 
characteristic for a certain material ||179|| . The high pressure limit of a cylindrical cell 



is set by the material strength and the ratio of outer-to-inner diameter ra/rj. In order 
to use the full size of a neutron beam the inner diameter has to be of the order of 
some centimetres. This easily results in large wall thicknesses leading not only to a 
strong background but also to an attenuation of the incoming beam and to rather 
difficult corrections. The background signal, however, not only depends on on the wall 
thickness, i. e. from the number of scatterers in the beam but also from the scattering 
cross-sections of the material. Corrections are difficult to perform if the attenuation of 
the beam is caused by the energy dependent absorption. Consequently, small ratios are 
used to keep these effects as small as possible. 

High pressure cells based on high strength Aluminum alloys have many of the 
required/desired properties. Vessels with small wall thicknesses withstanding high 
pressures produce a rather low background and nearly no absorption due to the low 
scattering and absorption cross-section. Bragg scattering of Al occurs only for wave 
numbers Q > 2.7 A~^. 

The experiments reported here were performed in a standard high pressure cell 
made of the Aluminium based alloy A17049.T6 of the ILL. The ratio of outer-to-inner 
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diameter of 2 allows to attain pressures up to 250 MPa at temperatures up to about 
340 K. The composition and scattering cross-sections of the high pressure material are 
listed in table |^. 

Table 3. Composition, scattering and absorption cross-sections of the high pressure 
material A17049.T6 of the vessel. 



Element 


wt% 


crc/[barn] 


<7i 1 [barn] 


(Tf/[barn] 


(Ta/ [barn] 


26 Fe 


0.3 


11.22 


0.40 


11.62 


2.56 


i4Si 


0.4 


2.16 


0.004 


3.49 


0.171 


24 Cr 


0.15 


1.66 


1.83 


4.35 


3.05 


30 Zn 


7.6 


4.05 


0.077 


4.131 


1.11 


12Mg 


2.5 


3.63 


0.08 


3.71 


0.063 


29 Cu 


1.5 


7.49 


0.55 


8.03 


3.78 


I3AI 


87.55 


1.49 


0.008 


1.503 


0.231 


A17049.T6 


100 


1.865 


0.026 


1.892 


0.358 



The transmission of the pressure cell was about 85 % and it was handled like a usual 
container. Restricting ourself to large scattering angles implying large wave numbers Q 
multiple scattering effects were neglected. The approximate correctness of this approach 
is demonstrated in figure where the correlation functions measured with and without 
pressure cell are compared. Apart from a tiny difference in the plateau height the curves 
are identical. 

Due to crystallization tendencies of OTP we had to use ^He gas as pressure 
transmitting medium. Any other attempt to use a liquid or OTP itself failed. With 
a scattering cross-section crcoh= 1-34 barn ^He gas is nearly transparent to neutrons. In 
addition the analysis is not complicated by the fact that a liquid pressure transmitting 
medium will show the same effects of dense liquid dynamics one is looking for in the 
sample. Clearly, '^He diffuses into OTP with a time constant of several hours. However, 
although the solubility increases with temperature and increasing pressure according to 
Henrys law, with a typical solubility of about 10~^ mol/ (gMPa) in aromates there is less 
than one He atom in the first coordination shell. 



1.2. Static Structure Factor 

The static structure factor of OTP-di^ was measured on the instrument D7, now used 
as a diffractometer without energy analysis. The pressure and temperature dependence 
of S{Q) is shown in figure S{Q) was obtained in absolute units by comparison with 
a vanadium standard. 

Under temperature variation at constant pressure (isobaric) the peaks behave as at 
ambient pressure. With decreasing temperature the scattering intensity in the Q —>■ 
region reduces due to the decreasing isothermal compressibility xt- S{Q = 0) = nxrkBT 
compares well with the values calculated from temperature and pressure dependent 
specific volumes [^. The peak positions shift slightly to higher Q values due to an 
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Figure 30. (*) Comparison of the correlation functions obtained without (open 
symbols) and in the high pressure cell (closed symbols) at Q = 1.2 (□,■), Q = 1.5 
(0,4) and Q = 1.8 A^^ (o,»)- The curves are nearly identical which validates the 
simplified approach of treating the high pressure cell like an ordinary container. 

increase in density n. Concerning the peak heights the maximum at 1.4 is not 
much affected, the maximum at 1.9 is enhanced and the prepeak is lowered when 
coohng. 

The effect of increasing pressure at constant temperature (isothermal) on S{Q) is 
similar to the one of decreasing temperature: S{Q — > 0) decreases and the peak positions 
move to higher Q. The intensity of the prepeak and the first maximum decreases while 
that of the second maximum increases. This behaviour is in agreement with results on 



other materials I [1801, 181 



Two other types of measurements are interesting in the light of MCT since the 
static structure factor is considered to be the trigger of the slowing down of the a- 
relaxation: isochoric measurements of S{Q), i. e. at (T, P) combinations leading to the 
same bulk density and isochronous measurements, i. e. at (T, P) combinations where 
the a-relaxation time is identical. Such measurements are included in figure |3l|. While 
along isobars and isotherms the static structure factors S{Q) evolve continuously, it is 
nearly identical along an isochor and an isochron. 

7.3. Square Root Singularity 

When the pressure dependence of the elastic scattering at constant temperature is 
monitored a behaviour slightly different as compared to the temperature dependent 



experiments is observed. In figure ^ left the isothermal DWF for incoherent scattering 
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Figure 31. (*) Static structure factor S{Q) under isobaric {P = 55MPa), isotliermal 
(T — 339K), isochoric {n = l.lOlgcni^'^) and isochronic (tq ~ 15nsec) conditions. 



(MoBbauer-Lamb factor) as a function of density is shown for several wave numbers Q. 
At high pressures or densities, in the glassy phase, the elastic intensity remains rather 
constant. Thus, no harmonic evolution had to be taken into account as in the case of the 
temperature dependent measurements. Releasing the pressure, an anomalous decrease 
of the elastic intensity sets in and indicates the presence of a faster relaxation. 

When constant pressure is applied the anomalous decrease shifts to higher 
temperatures with increasing pressure as apparent in figure right which shows the 
isobaric DWF as a function of temperature for three different pressures at Q = 1.8 A^^. 

Inspecting the isobaric decrease of fq in figure 32 at high temperatures more closely 
it seems that the slope of the decrease diminishes with increasing pressure and the 
transition range becomes broader. This slope is related to the strength of the fast 
f5 relaxation. One might speculate that the effect vanishes under extreme conditions. 
However, when the Q dependence of the elastic scattering is analyzed in terms of the 
mean square displacement (r^(T)) extracted from the Debye Waller factors in a Gaussian 
approximation (inset of figure [3^) this effect disappears. While around 200 K the (r^(T)) 
are almost identical the anomalous increase shifts to higher temperatures. The pressure 
dependence of the mean square displacement can even be scaled to a master curve taking 
the thermodynamic shift dTg/dP as the scaling parameter. Similar observations were 
made in the polymer polybutadiene ||183| , |184|| while in polyisoprene the application of 
pressure separates an additional dynamic process ||184| . 
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Figure 32. (a (*)) Plateau height fq of isothermal tagged particle correlations as 
a function of density at constant temperature T = 301 K for several Q-values. Solid 
lines are fits with a square root law leading to a critical density of ric ~ 1.100 gcm~^ 
which corresponds to a critical pressure of Pc — 42MPa. Towards smaller densites 
the integrated a intensities are independent of pressure. The inset demonstrates the 
square root law in linearized form, based on |182| . (b (*)) Incoherent elastic scattering 
intensity (M6f5bauer factor) as a function of temperature for different pressures at 
Q — 1.8 A^^. The linear behaviour at low temperatures is due to classical phonons 
exp(— Q^(r^(r))) with {r'^{T)) cx T. With increasing pressure the anomalous decrease 
shifts to higher temperatures. Above a temperature Tc{P) quasielastic broadening 
from structural relaxations comes into play and the integrated a intensity has to be 
taken into account. The inset shows the temperature dependence of the mean sqaure 
displacement which contains the complete Q information. 



7.4-- Critical Correlations 

We have already seen in the previous section that a kink in the Debye-Waller factor 
reveals the onset of the /3-relaxation. The anomalous increase of the DWF above a 
certain density indicates that the [3 process is also strongly influenced by the pressure 
or density. This fast process is responsible for the decrease of elastic intensity on 
costs of inelastic intensity. Again the /3-relaxation can be accessed directly on a TOF 
spectrometer like IN6. The spectral density below about 2meV varies strongly with 
pressure. In this region the dynamic structure factor satisfies the factorization rule for 



all pressures as already demonstrated in figure |T8. 

The /5-relaxation has been studied at variable pressure up 240 MPa | |182| , p.85|| . The 
tagged-particle correlations in the compressed liquid also decay in two steps and can be 
described by the MCT scaling laws. Unfortunately, the used instruments IN16 and IN6 
do not overlap (see table ^) and the dynamic range of IN6 alone is clearly not sufficient 
to warrant a four parameter fit. To get the pressure dependence of the /3-relaxation 
amplitude hq and the cross-over times to- we use the results from the a-relaxation 
analysis to be presented in the following section in more detail. Besides the DWF fq 
we fix the line shape parameter to the known value A = 0.77, which is compatible with 
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the frequency dependence u!^~"'{a = 0.3) shown in figure |T^. The DWF was determined 
from a Kohhausch fit to the long time tail of the low pressure data with fixed (3 = 0.6 
and tq{P). The insensitivity of line shape parameter A to pressure is supported by 
broad band light scattering measurements on the glass formers Cumen ||178|| and Salol 



186|| where in isothermal and isobaric experiments identical line shapes were found. 
Figure 133 shows isothermal pressure dependent tagged particle correlation functions 
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Figure 33. (*) Isothermal intermediate scattering function S{Q,t) at 316 K obtained 
on IN6 at three different wave numbers Q = 1.2, 1.5 and 1.8 A"^ in the /^-relaxation 
regime. Lines are fits with (^) with fixed line shape parmeter A = 0.77 and fixed 
plateau values (indicated by the thick circles). The pressures are from bottom to top: 
0,1 MPa, 15MPa, 45MPa, 60MPa, 90MPa, 150 MPa and 200 MPa. 

together with fits to (|l^). The correlation functions at the highest pressures decay to a 
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plateau which is constant in time and higher than the value for the liquid state: we are 
in the glassy state. For the lower pressures one recognizes a two step decay as in the 
temperature dependent case. While the first decay on the sub-picosecond time scale 
due to the microscopic dynamics is not covered by (|1^ it consistently describes the 
bending of ^{Q,t) into and out of the intermediate plateau fq over about one decade 
depending on the pressure. 

On approaching the cross-over from the liquid side, fg varies only weakly, whereas 
Hq and ta should become singular. In figure |3^ the power laws ([T8|) are tested for 
two isotherms by plotting Hq and )f:~^° vs. P, with a = 0.295 determined by A. 
Asymptotically, a linear behavior is observed allowing for a consistent determination 




P(MPa) P (MPa) 



Figure 34. /3-relaxation amplitude Hq (a) and cross-over time t^, (b) for two 
isotherms 316 and 324 K and three wave numbers Q = 1.2(n), 1.5(») and 1.8 A^^ 
{()). The data are hnearized as H q an d t^^°' according to (p^); hues are fits to Q 
averaged values. Reproduced from ]182| . 

of the mode-coupling crossover line Pc{T) or Tc{P). 

Together with T^. at ambient pressure and the P^ values obtained from the elastic 
scans six points of the dynamic phase boundary Pc{T) are obtained as summarized 

in table |. The (Pc,Tc) line has a slope dT^/dP ~ 0.26 0.28K/MPa which is in 

remarkable accord with the slope determined from a-relaxation, viscosity and Tg(P) 
data: the mode-coupling crossover line is parallel to lines of equal a response (see 
figure ^Ol). Lines of constant densities, however, have far steeper slopes of about 
0.7K/MPa. It confirms that the temperature influences the dynamics not only via 
the free volume. This finding is further illustrated by a direct comparison of isochoric 



correlation functions in figure The loss of correlations is the faster the higher the 
temperature. 

The above results imply that the dynamic crossover can be described by a single 
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Table 4. Critical pressure Pc and corresponding density Uc as function of temperature. 
The product Tc = na^(e/T)^^^ is found to be constant (its errors are estimated from 
ATc = ±3K and AP^ = ±5MPa). The choice of ct = 7.6 A and e = 600K is explained 
in the text. The values in the lower part stem from temperature dependent elastic 
scans where the critical temperatures can only be estimated. Nevertheless, the (T, P)- 
combinations are compatible with the phase boundery. 



T{K) 


Pc (MPa) 


ric (nm ■^) 




290 


0.1 


2.840 


1.495(4) 


301 


42(3) 


2.875 


1.499(6) 


316 


90(7) 


2.906 


1.497(5) 
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123(7) 


2.931 


1.501(9) 
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Figure 35. (*) Isochoric tagged particle correlation function $(Q, t) for Q=1.8 A ^ at 
a density n ~ 1.083g/cm'^. The {T, P) combinations are from top to bottom: (292,0.1); 
(300,12); (310,26.5); (319,39) in K and MPa. Lines are fits with a fixed t^(r). The inset 
shows the results for the /3-relaxation amplitude Hq in linearised form for three wave 
numbers all extrapolating to a consistent single value of rtT^^/'* cx F. For explanation 
of F see text. 

separation parameter a{T, P) which, however, does not depend on density alone. 
Instead, it appears that density n and temperature T can be condensed into an effective 
couphng r oc nT"^/^ which is constant within experimental error as shown in table ^. 
Furthermore, constraint fits to all individual ^{Q, t; T, P) confirm that the temperature 
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and pressure dependence of the amplitude and the cross-over time can be described by 
the power law (|T8|) in conjunction with just one separation parameter a = (T — Tc)/Tc 
(figure 35). At least in a certain neigbourhood of Pc{T) this parameter is linear in the 
effective coupling T oc nT~^/^ which is known to characterise equilibrium properties 
of a dense soft sphere fiuid with repulsive core potentials of the Lennard- Jones 
type ||187|| . Note, that a plot log^c vs. logTc yields a slope of 0.27 very close to the 
theoretically expected value 1/4 for k = 12. 

In a supercooled and compressed liquid, we are actually in a high-density regime 
where the static structure factor S{Q) is determined mainly by the repulsive part of 
the potential ||187| , |188i |189|| . Since this repulsive part can be modeled well by a soft 



sphere system our results are thus evidence that for the glass transition subtle changes 
in the static structure factor are responsible for the slowing down of the dynamics. This 
view is also supported by computer simulations of binary soft sphere ||190| , |191| , |192 



and Lennard- Jones ||148|| systems, in which the glass transition was found to occur at 
constant Tc- It is astonishing that these results apply literally in a complex molecular 
liquid like OTP. In fact, the comparison with a Lennard- Jones liquid can be made 

almost quantitative. To do so, we need the two parameters of the 6 12 potential: 

a van der Waals diameter a and a potential depth e. cr = 7.6 A is estimated from 



structural information [Q. This is in good accord with the prepeak in the static 
structure factor at Q ^ 0.85 . . . 0.9 A~^ ^2], |157|| and the scaling factor needed to describe 
the Q dependence of the critical Debye-Waller factors fg and the amplitudes hg in 
figure 0. A potential depth e = 600 K has been determined by calibrating a three-site 
OTP model to experimental density and diffusivity |135|| . Using this input, we obtain 
Fc = 1.50 (table ^) which is closer to the Lennard-Jones result F^'^ ~ 1.25 |[146| , |148 
than one could reasonably expect. 

Bengtzelius | 146 | calculated Tc{n) for different densities for a Lennard-Jones 
system. In order to compare our Tc{n) with his result we rescale his axis such that they 
coincide at Pq. Our slope (5.02 10"'^ gcm~^) is in good accord with the theoretical one 
(6.22 10~^ gcm~^). From the rescaling we get the two Lennard-Jones parameters: the 
potential depth e/ks — 608 K which astonishingly is the value found in the simulations 
||135|| and the diameter cr ~ 7A being in good accord with the van der Waals radius. 

We have found in the pressure and temperature dependent measurements of the 
static structure factor that along an isochronous line (r^ = const.) S{Q) does not 
change significantly (figure 0). According to the theory the spectra S{Q,lj) should all 
be identical (apart from the possibility of a trivial shift in the overall time scale to)- Iii 
order to test this behaviour directly we measured the dynamics for several combinations 
of temperature and pressure leading to the same F (and not as usual along isothermal, 
isochoric or isobaric paths). A direct comparision of the dynamic structure factors in 
figure ^ shows indeed that the sprectra are identical within statistical accuracy. This 
implies that the plateau height fg, the /3-relaxation amplitude hq and the cross-over 
time tcr are identical along a constant F, i. e. along a constant Xa-line. Note, that for 
our purposes it is sufficient to compare S{29, uj) after subtracting the empty can spectra. 
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Figure 36. (*) Incoherent dynamic structure factor S{29,uj) at a fixed scattering 
angle 26 = 123° for several combinations of temperature and pressure leading to the 
same value for F. Note, the spectra are identical within experimental error. 



Table 5. Connection between expansivity a, compressibility (3 and the repulsive part 
of the potential expressed via (^) along the dynamic phase boundery i. e.. 



T(K) 


Pc (MPa) 


alOV K-i 


PIO^/ MPa-i 




290 


0.1 


7.2 


5.0 


0.31 


301 


42(3) 


6.4 


4.4 


0.28 


306 


60(3) 


6.1 


4.1 


0.26 


316 


90(7) 


5.7 


3.8 


0.25 


317 


100(5) 


5.6 


3.7 


0.24 


324 


123(7) 


5.3 


3.6 


0.24 



Interestingly, a relation between the effective coupling T or, more exactly, the 
repulsive part of the potential and some macroscopic quantities can be established. 
To do so we assume the specific volume v{T, P) to be linear dependent on temperature 
and pressure via the thermal expansivity a and the isothermal compressibility (3: 

v{T, P) = v{T,, 0)(1 + a{T - TJ - (3P,{T)) (28) 

with Tc the critical temperature and Pc(T) the critical pressure. At the phase 
boundery = fT^/^ = const., thus we can derive the following equation using 
T = Tc(l + (T — Tc)/Tc) and developping in a Taylor series up to first order 
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where we have used the fact that the phase boundery follows 

Using our (Tc, Pc) combinations, the values for the expansivity and compressibility at 
these @ and dT/dP\r ^ 0.28K/MPa gives us the result in table ||. Astonishingly, 
equation (^Of ) holds quite good indicating a deeper link between the microscopic repulsive 
part of the potential and the macroscopic quantities expansivity and compressibility. 

It would indeed be interesting to check on other fragile glass froming liquids, 
whether the critical mode-coupling phase boundery ist parallel to lines of constant 
a-relaxation time and whether an effective coupling and a link to the expansivity and 
the compressibility can be established in a similar manner as we found for OTP. 

7.5. a -Relaxation 

On the BS instrument IN16 we measured the slowing down of the structural a-relaxation 
as a function of pressure. In a first survey, we could show that the mean relaxation 
time follows the pressure dependence of the viscosity. The shape of the a-relaxation 
spectra was then determined in some long-running scans. At three different points 
(T, P) of equal viscosity, the quasielastic broadening of all spectra coincide as shown in 
figure ^ for three values of Q. The microscopic relaxation time is proportional to the 
macroscopic viscosity rj and the spectral distribution is quite independent of pressure. 
This is confirmed from photon correlation and dielectric measurements ||6^ . 



By construction of a master curve for isothermal spectra at 320 K we can further 
illustrate the invariance of the line shape and the viscosity scaling. Motivated by the 
result of figure ^ we use viscosity data ||2^ to rescale the times t — > t = tri{Po) /rj{P) with 



Pq = 0.1 MPa similar to the procedure described in section |6.5| . The rescaled 
coincide over two decades as shown in figure ^ for three Q values. The scaling is well 
satisfied. The a-relaxation amplitude and the line shape parameter f5 as derived from 
the iterative procedure are shown in figure Both parameters are nearly independent 
of pressure, which in turn justifies the analysis of the /3-relaxation presented in the 
previous section. The strechting parameter /3 is similar to that derived from previous 
investigations at atmospheric pressure |^ : The distribution of relaxation times is pretty 



much uneffected by the pressure and temperature changes. The mean relaxation times 
(r(P)) as obtained from Kohlrausch fits with fixed {Aq)p and jS and after elimination of 



the Q dependence are also shown in figure ^ It follows the same pressure dependence 



as the viscosity, in agreement with the isoviscous measurements. 

Thus, a time-temperature-pressure superposition principle holds. It is well 
established by PCS and dielectric spectroscopy in OTP as well as in some other fragile 
liquids ||70|, |193|| . Our measurements extend the validity of this scaling principle down 
to the picosecond time range. 

The pressure dependence of isochronous relaxation times is compiled in figure 
Lines of isochronous a-relaxation times are connected by the same slope of about 
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{dT/dP)r — (0.28±0.01) K/MPa over about twelve orders of magnitude regardless of the 
time scale. At long times the pressure dependence of the glass transition temperature 
(dTg/dP) is the same as in the second to microsecond time range covered by 

photon-correlation spectroscopy |^ and dielectric spectroscopy down to the sub- 
nanosecond regime covered here by neutron scattering. 




tico (|J.eV) 



Figure 37. (*) a-relaxation spectra from IN16 at Q = 1.2,1.5 and 1.8 A^^ for 
three combinations of temperature and pressure leading to approximately the same 
viscosity. The spectra are scaled to their value at w = 0. The hached area represents 
the measured resolution function. The {T,P) combinations are: 306K, 0.1 MPa (□); 
320 K, 50 MPa (■) and 335 K, 100 MPa (0). 
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Figure 38. (*) Master curves $q(£)/j4q(P) for the isothermal a-relaxation of the 
tagged particle fluctuations at 320 K and different pressures, obtained by rescaling 
times with viscosity r]{P). 
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Figure 39. (*) Kohlrausch parameters /3q (a), Aq (b) and (r) (c) as obtained from the 
scaling analysis as a function of pressure. In each fit, two parameters were fixed at the 
values obtained from the master curve analysis. Amplitude Aq and stretching exponent 
/3q are roughly independent of pressure in the liquid state. The mean relaxation time 
(r) scales with visosity (line). 
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Figure 40. (*) Pressure and temperature dependence of isochronous a-relaxation 
times. Note that from Tg corresponding to relaxation time of the order of 100 sec 
down to nsec all slopes are identical. In particular the dynamic phase boundary Tc{P) 
is parallel to lines of equal a response. Constructed using values from |4^, |60[ |. 
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8. Vibrational Dynamics [36] 

When the dynamic anomahes of the glass transition are discussed it is particularly 
important for molecular glasses, which have a large number of internal degrees of 
freedom, to study at least qualitatively the vibrational mode structure. This seems 
to be even more important in view of the fact that low-frequency vibrational features 
- for obscure reasons called Boson peak - persist well into the supercooled phase. 

This can be done by comparison of the absolute frequency scales of the vibrational 
density of states between the ordered and disordered state. From single crystal 
dispersion measurements - if available - information about low-energy internal modes 
and hybridisation of external and internal modes can be obtained. 



8.1. Single Crystal: Phonon Dispersions 

The phonon dispersions of single crystal OTP at 200 K along the three main symmetry 



directions [100], [010], and [001] are presented in figure |4l|. Some optic-like phonons are 
detected as well, as an example, the lowest optic branch in [010] direction is included 



in figure ^ However, for a detailed study of them thermal neutrons instead of cold 
ones would be needed. In the small Q limit the sound velocities of acoustic excitations 
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Figure 41. Phonon dispersion relations of crystalline OTP as measured on the triple 
axis spectrometer IN12 at 200 K. The wave vectors are given in reciprocal lattice units, 
excitation energies in THz (ITHz = 33.36 cm~^ = 4.1357meV). Taken from p^ . 

obtained from the initial slope are summarized in table |^. In all lattice directions, 
longitudinal sound modes are almost twice as fast as transverse modes. The single 
crystal data may also be compared with sound velocities in the glass. Results from 
Brillouin light scattering are summarized in table 0. Taking the simple arithmetic 
average over the three crystal axes, the mean longitudinal and transverse sound velocities 
(■^l.t) exceed those of the glass by about 20 % and 27 %, respectively. 
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Table 6. Sound velocities (in km/s) for three phonon branches and different lattice 
directions, obtained from linear fits to the low-Q limit of the measured phonon 
dispersions. 

"[100] «[010] W[ooi] (v) {v~^) 

VL 3.71 2.68 3.30 3.23 3.11 
VT, 1.82 1.52 1.97 

1.75 1.71 

VT2 1-67 - - 



Table 7. Sound velocities in the glass as measured by Brillouin scattering, using 
visible light or X-rays. 
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2.70 
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light 
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1 


117 
119 
121 





8.2. Glass versus Polycrystal: Density of Vibrational States 

Without crystalline order, it is no longer possible to measure selected phonon modes 
with well-defined polarization and propagation vector. For polycrystalline or amorphous 
samples, the distribution of vibrational modes can be conveyed only in form of a spectral 
density of states (DOS). 

The DOS of glassy and crystalline OTP at 100 K are shown in figure They were 
calculated directly from the spectra at the highest scattering angles S{26, u) without 
interpolation to constant Q. Multiphonon contributions were calculated be repeated 
convolution of g{uj) with itself and subtraction from S{26,uj) in an interative procedure 
described in detail in ||107|| . Comparing our results with model calculations ||194|| , we 
assign all vibrations below z/g to the 16 degrees of freedom needed to describe the crystal 
structure. For 16 low-lying modes in a molecule with 32 atoms, we expect an integrated 
DOS 

= 16 

dug{u) = - = 0.5. (31) 

Rather broad distributions are found for both phases showing a pronounced gap 
above z/g ~ 5THz (ITHz = 33.36 cm~^ = 4. 1357 me V). In the glass a first shoulder 
around 1.5 THz is followed by a second at 3.5 THz in accordance with results from Raman 



studies [115, 116|. As expected, the crystal DOS is more structured, in particular in the 



low energy region. Distinct peaks at 0.6, 0.8, 1.1 and 1.5 THz become apparent. They 



are due to strong contributions from zone-boundary modes (figure ^). Compared to 
the glass, significant density is missing in the low energy region and in the range from 1.5 
to 3 THz, and reappears at higher frequencies around 3.5 THz. The gap for z/ > 5 THz 
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Figure 42. (a) Density of vibrational states of glassy and crystalline OTP at 100 K. 
(b) To emphasize the excess density of states of the glass over the crystal g{v)/v'^ is 
shown. The arrows mark the Debye limit — > calculated from sound velocities and 
densities at 220 K for the glass and at 200 K for the crystal. Taken from js^ . 



common for both states and indicates a separation from high energy intramolecular 
modes. 

The DOS of the polycrystalline sample is in accord with the measured dispersion of 
the single crystal: A small shoulder at 0.4 THz can be attributed to the transverse 
acoustic zone-boundary phonons in [100] direction, and the main peak at 0.6 THz 
corresponds to the transverse acoustic zone-boundary phonons in the other two lattice 
directions. The peak at 0.8 THz reflects the longitudinal acoustic zone-boundary phonon 
in [001] direction and the transverse optic in [010] direction. 



In order to show the low energy modes on an enhanced scale, we plot in figure 42 



g(v)jv^ which in the one-phonon approximation is proportional to S{Q^v) itself. In 
this representation, the excess of the glass over the crystal becomes evident. A well 
defined frequency peak appears around 0.35 THz which is down shifted with respect to 
the first peak of the crystal at 0.6 THz and superposed to a long tail which is similar 
for both, glass and crystal. Note that the maximum of the boson peak at 0.35 THz is 
located below the lowest acoustic zone-boundary phonons in the crystal. 

Another dynamic observable which can be obtained from neutron scattering is the 
atomic mean square displacement (r^(T)). Roughly speaking, (r^(r))'^/^ measures the 
volume to which an atom remains confined in the limit t ^ oo. For a large class of 
model situations (harmonic solid, Markovian diffusion, . . . ) it can be obtained directly 
from the Gaussian Q dependence of the elastic scattering intensity 

^(Q,z. = 0) = exp(-g2(r(T)2)). (32) 

For a harmonic solid. 
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(with f3 = hu/k^T) crosses over from zero-point oscillations (^'^(O)) to a linear regime 
{r^{T)) oc T. Note, for a meaningful application of (|33|), we have to assume the rigid- 



body limit (Sect. 3.1, Ref. [107]). Therefore, M has to be taken as the average atomic 
mass which in the case of OTP is 7.2. In any real experiment, which integrates over the 
elastic line with a resolution Az/, one actually measures atomic displacements within 
finite times t^ ~ 27r/Az/. 



Figure ^ shows mean square displacements of OTP, determined according to (p^) 
from elastic back-scattering (fixed window scans on IN13, with t/^ ~ 100 — 200 psec) 
and from Fourier-deconvoluted time-of-flight spectra (taking the plateau S{Q,t^) with 
— 5 — 10 psec from Fourier transformed IN6 data). For the glassy sample, a direct 
comparison can be made and shows good agreement between IN6 and IN13. 
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Figure 43. Mean square displacement (r^(r)) of glassy OTP from elastic back- 
scattering (IN13, 0) and from time-of-flight spectra (IN6, plateau value of S{Q, t), □). 
For comparison, IN6 data of polycrystalline OTP are also shown (■). The lines are 
calculated from the densities of states at 100 K of the glass (dashed) and the polycrystal 
(solid line). The inset shows the relative value of d — (r^)^/^ at 100 K for the crystal 
and the glass when the integral (^) is restricted to modes with < ly < i^' . Taken 
from Hi. 

The lines in figure ^ are calculated through (|33D from the DOS at 100 K. For 
low temperatures, the (r^(T)) are in full accord with the values determined through 
(P^). This comparison can be seen as a cross-check between the analysis of elastic and 
inelastic neutron scattering data. 
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Equation ( p3D not only gives the absolute value of (r^(T)), but enables us also to 
read off which modes contribute most to the atomic displacement. To this end, we 



restrict the integration (|33D to modes with < u < u'. The inset in figure i3 shows 
the relative value {r^{iy';T)y/^ / {r^{iy;T)y/^ for T = 100 K as function of z/'. Modes 
below 0.6 THz in the crystal, or 0.4 THz in the glass contribute about 55% to the total 
displacement; 90% are reached only at about 2 THz. This means that the modes which 
are responsible for the mean square displacement and the Debye- Waller factor are not 
rigid body motions alone, but contain a significant contribution from other modes. 

Anyway, one immediately recognizes that for all temperatures the mean square 
displacement in the glass is larger than in the crystal. For higher temperatures the 
measured mean square displacements are clearly larger than expected from harmonic 
theory. These anharmonic contributions to (r^(T)) start around 140-150 K. Around the 
same temperature deviations from the proportionality (r^(T)) cx T can be observed in 
coherent elastic scans and in the MoBbauer experiment ||162| (figure |15| in section |6.3|) . 
In addition, around 140-150 K the temperature dependence of the experimental specific 
heat cp(T) |^ shows a slight change for both glass and crystal. 

Anharmonicities can directly be observed via the temperature dependence of the 
density of states. With increasing temperature the first moment, which provides a 
measure for the center of gravity, clearly shifts to lower frequencies for the glass and 
even more for the crystal. Note however, the strong increase above 240 K in the glass is 
connected with the glass transition. 

Our vibrational density of states can also be checked againts the heat capacity, 
which for a harmonic solid is given by the integral 

'sinh2(/3/2) 

With a Debye DOS, this yields the well-known cp oc T^. Therefore, in figure ^ 
experimental data for the specific heat of glassy and polycrystalline OTP are plotted 
as cp/T^. In this representation, a boson peak at 0.35 THz is expected to lead to a 
maximum at about 4K. 

The lines are calculated through (|3^) from the neutron DOS. They agree with 
the measured data for both glass and crystal in absolute units and over a broad 
temperature range. Similar accord has been reported for a number of other systems 
|TT5i |T95i |T96| , |T0|, |197 . 



Cp(T) ~ cv(T) = N^tR I dz/ (7(^^)-^7T7777;^ • (34) 



At higher temperatures, the heat capacities of crystalline and glassy OTP differ 
only little {e.g. at 200 K: Cp^'* = 182.8 Jmor^ K~i and cf.''''' = 186.1Jmol-i K'^ [|3[). 
Towards higher frequencies, the DOS becomes less sensitive to the presence or absence 
of crystalline order and remaining differences (clearly visible in figure ^2|) are largely 
averaged out by the integral (pi]). 

Altogether, these observations suggest that the mode structure in the glass and the 
crystal is quite similar except for an overall softening in the glass. Concerning the above 
presented analysis of quasielastic scattering in the supercooled hquid we can state that 
above 0.6 THz a very large portion of inelastic intensity in S{Q, v) of glassy OTP must 
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Figure 44. Experimental heat capacities Cp/T^ of glassy (□) and crystalline (■) 
OTP compared to cp/T^ calculated from the density of vibrational states g{i') at 
100 K of the glass (dashed line) and the polycrystal (solid line). Taken from ||36|| . 



arise from optic modes. Hybridisation and coupling between inter- and intramolecular 
modes plays an important role for frequencies higher than 0.6 THz [|198| , |115|| . As a 
consequence, the quasielastic scattering, which is confined to below 0.25 THz (ImeV), 
is clearly dominated by rigid-body motions and the analysis in terms of relaxations 
remains reasonable. 
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9. Conclusions 

As the reader has seen in the preceeding sections, neutron scattering experiments on 
the fragile van der Waals liquid Orthoterphenyl provide clear evidence for a dynamic 
cross-over from a liquid-like to a solid-like motional mechanism occuring around the 
critical temperature Tc. 

Mode-coupling theory gives a quantitative consistent description of the dynamics in 
the supercooled state. Time-temperature-pressure superposition, known from Hz-MHz 
spectroscopy, is valid down to the picosecond time range, where macroscopic properties 
of the material evolve from their microscopic dynamics. Spectral line shapes observed by 
coherent and incoherent neutron scattering experiments are invariant under temperature 
and pressure variations. The mean relaxation time is proportional to the viscosity and 
states of equal kinetics are connected by the same slope dT/dP, regardless of the time 
scale. 

Temperature and pressure dependent Debye-Waller factors show an anomalous 
decrease that can be described by a square root law. Whether or not the square root 
law results in a sharp singularity is however sometimes not compelling, in particular 
for small wave numbers Q or at the structure factor maximum, i. e. where the effect 
is expected to be very small anyway. The lost scattering intensity reappears as an 
quasielastic process and is identified with the /3-relaxation whose time correlations 
follow the predicted scaling behaviour. More importantly, all correlation functions can 
be described by the same MCT (3 correlator ( p!7D with one consistent set of exponents. 
The fit parameters show the predicted Q dependence and consistently converge to the 
same critical Tc{P). They allow the unambigous determination of a dynamic crossover 
line Tc{P). This line is parallel to the lines of equal viscosity: a- and /3-relaxation 
are driven by the same parameter. At least in a certain neighbourhood of Tc{P) this 
parameter is linear in the effective coupling T oc nT~^/^ which in turn ultimately depends 
on the short range interaction of the molecules. Our experiments have also allowed for 
to separate temperature from density effects and to relate at least qualitatively dynamic 
and structural changes. 

Furthermore, we have excluded any intramolecular origin of the observed anomalies 
by measurement on partially deuterated samples and by comparison of frequency 
scales of the glass and the crystal. Any structural reason is excluded by the smooth 
temperature and pressure evolution of the static structure factor S{Q). Also in that 
sense OTP represented a good model system for testing the predictions of MCT. 

However, mode-coupling theory does not describes some other important aspects 
like slow processes in the vicinity of Tg, slow secondary Johari-Goldstein /^-relaxation 
5^, cooperativity, dynamic heterogeneity | [L99| | and the vibrational dynamics. We 



should also mention that in a molecular liquid the first critical exponent a has never 
been observed purely and that the appearance of the so-called " knee" below the critical 
temperature is still an open question. 

Experimentally, the identification of the dynamic transition is far from being trivial 
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and the comparison of neutron scattering and MCT predictions still remains qualitative 
in many aspects. Complicated spectral shapes have to be measured with high precision, 
a good Q resolution and over a large time or frequency window. We presented two 
different approaches that partly satisfy these demands: the combination of different 
spectrometers and the analysis using mastercurves. 

On a fundamental level, the presence of a de Gennes type narrowing in 0TP-di4 
reinforces our very basic assumption that incoherent and coherent neutron scattering 
probe self and density correlations of OTP molecules, respectively. Of course, we 
measure correlations of single H, D and C atoms which, in addition to a translational 
component, also contain contributions from rotational motion relative to the center-of- 
mass. Whether such contributions are relevant and more sophisticated MCT extentions 
have to be used can only be decided by joint use of different experimental techniques 
including molecular dynamics simulations with appropriate models. 

For the time being we believe that in our study we have exploited the very limits of 
neutron scattering in its present state of the art. We have also demonstrated that some 
very interesting experiments like the separation of coherent and incoherent scattering 
are feasible, but better and more accurate results have to await new high-flux neutron 
sources. The total efficiency has to be increased by at least one order of magnitude 
by increasing the incomming flux and / or by a much higher spectrometer efficiency. 
Upgrading and development of all types of spectrometers is currently under way. 

Other methods, e. g. synchrotron radiation, optic Kerr effect, impulsive stimulated 
light scattering or high frequency dielectric spectroscopy, among others, seem also to be 
very promising for the investigation of glassy dynamics in particular as they couple to 
different degrees of freedom. They may open the way to a comparison of experimental 
data with more sophisticated versions of the theory. 
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